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Abstract: It is widely believed that an atom interacting with the electromagnetic field (with total initial energy 
^ well-below the ionization threshold) relaxes to its ground state while its excess energy is emitted as radiation. 

f~i . Hence, for large times, the state of the atom+field system should consist of the atom in its ground state, and 
Q-l a few free photons that travel off to spatial infinity. Mathematically, this picture is captured by the notion 
(— I ■ of asymptotic completeness. Despite some recent progress on the spectral theory of such systems, a proof of 
■ relaxation to the ground state and asymptotic completeness was/is still missing, except in some special cases 

a (massive photons, small perturbations of harmonic potentials). In this paper, we partially fill this gap by proving 
, relaxation to an invariant state in the case where the atom is modelled by a finite-level system. If the coupling 

to the field is sufficiently infrared-regular so that the coupled system admits a ground state, then this invariant 
' state necessarily corresponds to the ground state. Assuming slightly more infrared regularity, we show that the 
K*" , number of emitted photons remains bounded in time. We hope that these results bring a proof of asymptotic 
completeness within reach. 

vn- 
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1 Model and result 
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1.1 Introduction 



This paper fits into a broader project of rigorously controlling interacting Hamiltonian systems that exhibit 
irreversible behavior. From a more concrete point of view, the problem treated here is inspired by nonrelativistic 
QED which, in the last two decades, has proven to be a fruitful testing ground for mathematical techniques in 
d I quantum field theory (for an overview, see the book Il29ll ) Sacrificing precision for the time being, the setup is 
as follows. We consider a model of an atom interacting with the electromagnetic field (or, a scalar field, as we 
will assume in this paper for the sake of simplicity). The field is described by the Hamiltonian Hp on Hilbert 
space J^, it is non-interacting and describes freely propagating scalar bosons with a linear dispersion law (the 
polarization of the photons does not play a crucial role in the physical problem, hence we omit it from our 
model). The joint dynamical system (atom + field) is described by a Hamiltonian H on a Hilbert space J/f and 
the field is coupled in a nontrivial way to the atom. Furthermore, we assume that the atom cannot be ionized, 
i.e. in the absence of coupling its spectrum is discrete. We wish to address the long time behavior of the system. 
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1.1.1 Approach to a stationary state 

Let O be a local observable in (i.e. bounded operators on ,y/f) that should be thought of as probing the 

atom and the field in its spatial vicinity. Let 4*0 G ||^o|| = 1 be the initial state vector (0 refers to time 
i = 0) and '■— e^'*^^o is the time-evolved state vector. We consider the expectation value ( (•, •) is the scalar 
product) 

{0)t:^{^u0^t) (1.1) 

As explained in the abstract, we substantiate the claim that this expression converges to an asymptotic value, as 
< — > oo. i.e. {0)t — > (0)oo and that (0)oo is independent of If the Hamiltonian H admits a groimd state *gs, 
i.e. Egs = inf cr{H) is an eigenvalue, then this implies that (0)oo = (^gs, O^gs) since one can choose = ^gs- 
In particular, this means that the eigenvalue Eg^ is simple and that H cannot have any other eigenvalues. It 
is this last claim that has been established up to now in great generality. In fact, one can even prove that H 
has absolutely continuous spectrum, apart from the simple eigenvalue Egg, see (U |l2l [l3l [5]. However, this 
seems not sufficient to prove that dl.lb converges as t ^ oo, because e~'*^ appears twice in this expression 
(!). Let us immediately add that we find it not at all inconceivable that there is some easy way around this 
problem, allowing to apply the techniques used for the spectral analysis of H to determine the asymptotics of 
(|l.lb . However, this is not the strategy of the present paper. We prove the convergence to an asymptotic value, 
{0)t —?' (O) oo not by spectral considerations, but by exhibiting explicitly the irreversible density matrix evolution 
1^*0) (^o| If H is too infrared-singular to admit a ground state, then the convergence {0)t — > (0)oo 

can still hold provided that the observable O 'does not see' the low-energetic photons. In that case the asymptotic 
value (0)00 is a state (a positive, normalized functional in O) that is not of the form {0)oo = (*gs, O^gs). From 
the point of view of our technique, this case is no different from the infrared-regular case. 

Finally, we mention that the problem of 'return to equilibrium' at positive temperature has been studied 
with much more success since in that case, the problem can indeed be reduced to the study of the spectrum of 
an operator - the so-called standard LiouviUian - acting on an appropriate Hilbert space, see [7| for references. For 
our technique, there is no difference between zero and positive temperature, and the present result on approach 
to a stationary state was, up to some irrelevant details, in fact already contained in [7J. 

1.1.2 Scattering theory 

Let us again assume that H admits a ground state ^tgs (and no other eigenstates). The intuition is that the 
evolved state vector ^I'f should, at large times t, look like the ground state with a few free photons that travel off 
to infinity. 

To make this intuition more precise, one introduces an identification operator Jid : J^, such that 

Jid^i = 5'gs with n the field vacuum and in general Jjd maps photon state vectors 'ifp into state vectors that 
we think of as 'ground state together with a free photon wavepacket '^p'. One defines the asymptotic wave 
operators W^+, W~ by 

VF=^*F := lim e'*^^-^'^^') Jide-'*^''^'p' (1.2) 

for in a dense domain. Asymptotic completeness (AC) asserts that the operators : Jfp ^ ^ exist 
and are unitary. The main unproven aspect of this statement is RanW^^ = Indeed, what can happen in the 
presence of (massless) photons is that — e^'*^5'o will contain ever more (as t 00) photons with ever smaller 
energies. But if this is the case then = e'*^4't can clearly not equal VF+^'f for any £ -y^- For this reason 
it is important to bound the number of photons of the state vector . Indeed, if one eliminates the soft photons, 
either by making them massive or by introducing a sharp irtfrared cutoff in the coupling, then all questions can 
be answered and in particular AC follows, see [11. 16|. Our result provides a strong exponential bound of the 
form (^Pf , c'^^^E't) < C with k a small positive constant and C independent of time. We believe that with this 
bound as an input, a proof of asymptotic completeness is within reach (see e.g. iTTSl ) and we hope to pursue 
this in a subsequent paper. Apart from the case of massive bosons, asymptotic completeness can be established 
if the particle is a harmonic oscillator and the coupling to the field is linear, so that the full Hamiltonian can 
be explicitly diagonalized jj). Small perturbations (due to small anharmonicities in the particle potential) can 
be handled by an expansion introduced in |23|, see |28|. Of course, if one considers models not described by 
quantum field theory, for example; iV -particle quantum mechanics, more results are available and we refer to 
IITTl for more references and backgroimd. 
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1.2 Setup 



Let be a Hilbert space (modeling the small system) with a self -adjoint Hamiltonian Hg. The field is given by 
the one-particle dispersion relation \q\ and the Hamiltonian of the whole field is given by 

Hp := / dq\q\agaq (1.3) 

acting on the bosonic (symmetric) Fock space = r(t)) with t) — L^(K'*) the one-particle space. Here a*,aq 
are the creation/annihilation operators (actually, operator-valued distributions) of a mode with momentum 
q G M."^ satisfying the canonical commutation relations [aq,a*,] — d{q — q'). We refer to e.g. [9) for a review 
of these notions and precise definitions. The Hilbert space of the total system consisting of small system and 
field, is ^ = Ms ® J^, and we simply write iJs and for the operators Ha,®l and 1 ® Hp acting on . 
For simplicity, the coupling between field and the small system is assumed to be linear in the creation and 
annihilation operators and of the form XHi, with A e R a (small) coupling constant and 

Hi = D(g> f dq U{q)a* + , (1.4) 

for some G t) and Hermitian matrix D = D* E ^(J^). Our technique works equally well if one considers a 
finite sum of such interaction terms, or if one adds a sufficiently small quadratic interaction term, but we prefer 
to keep the setup as elementary as possible. The formal total Hamiltonian of the system is hence 

H := Hs + Hp + XHi, on Jf". (1.5) 

To construct H rigorously, we assume throughout that 

(0,(1 + ^)0),, - / dq(l + ^)|0(q)|2 <oo. (1.6) 

\<l\ Jr" \q\ 

which yields, by a standard estimate: 

||i?i*ir < 2||£'||(</',(l + l/k|)</')[, («',i?F*> < C(*,ffF^'), for^- e Dom(ffF) (1-7) 

It follows that Hi is an infinitesimal perturbation of Hp and by the Kato-Rellich theorem, the Hamiltonian H is 
self-adjoint on Dom(_ffs + ^f)- 

The following assumption is the key ingredient of our analysis, as it expresses that correlations of the free 
field vanish in time sufficiently fast. 

Assumption A (Decay of correlation functions). 

[ dt {1 + t)°'\h{t)\ <oo, withh{t):^[ dge-'*l''l|(/)(g)p, (1.8) 

for some a > 0. 

This assumption implies some infrared regularity of the model. In particular, if it is satisfied with a> 1, then 
the Hamiltonian H has a ground stat^ as one can establish by, for example, the techniques in IT7l[T8l l4l[T4l[T]|. 
In SectionO we give a condition on the form factor that is sufficient for Assumption lAlto hold. 

The next assumption is meant to exclude situations in which the atom is poorly coupled to the field. In 
particular, if the function <f>{q) vanishes identically, then one cannot expect relaxation to the ground state, and 
this should surely be excluded. We assume that a{Hs), the spectrum of Hs, is nondegenerate, and let Pe for 
e e cr{Hs) be the corresponding (one-dimensional) spectral projectors. We introduce the nonnegative numbers 

j-(e, e') := i:i[P^DP^,DP^] h{e - e'), (1.9) 

■'in fact, weaker conditions suffice for the existence of a ground state. For example, discontinuities in the 
form factor 4> can invalidate our results (see Section ll.S.lb , but the existence of the ground state depends solely 
on the behavior of (j) near 0. 
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where h{e) ~ e'*'^ft,(t) is well-defined by virtue of Assumption lAl It can be written in a more intuitive 

fashion as In J^^ dqS{\q\ — (e — e')) which one recognizes as the textbook Fermi Golden Rule expression 

for a scattering rate. One deduces that j(e, e') = whenever e' > e. Physically, this expresses that the field is in 
the vacuum state and it can only absorb (and not emit) energy. 

Assumption B (Fermi Golden Rule). We assume that the spectrum of is non-degenerate (all eigenvalues are simple) 
and we let eo := min(T(i?g) (atomic ground state energy). Most importantly, we assume that for any eigenvalue e e 
a{Hs),e ^ eo, there is a sequence e(i), i = 1, . . . ,n of eigenvalues such that 

e = e(l) > e(2) > . . . > e(ri) = eo, and Vi = 1, . . . , n - 1 : j(e(i), e(i + 1)) > (1.10) 

zvith as defined above. 

The numbers j(e,e') should be viewed as 'jump rates': We define the one-dimensional spectral projector 
Po = |^^)(^^| G with f2 the vacuum vector in the Fock space Mp. If the joint atom-field system is 

described by the density matrix po — ® Pn at time < = 0, then formal perturbation theory (Fermi Golden 
Rule) suggests that the probability to find the atom in state e' 7^ e at a later time i > 0, is 

Tr[Pe'e-"^Poe"^] = j(e, e'){\H) + 0{{\Hf) (1.11) 

The rigorous version of this formula, given in Proposition 13 .31 is a crucial ingredient of our analysis. 

1.3 Initial states and observables 

We now define the class of initial states po and observables O that we consider. For ^ G f), let yV(V') G ^(J^) be 
the Weyl operator 

W(V') = e'*W, $(^) := j dq {i,{q)a*+^a,) (1.12) 

We use K , >^ as labels to denote objects characterizing the initial state ('left boundary') and observable ('right 
boundary'). We pick '0x1 € f), an atom observable Os G ^{J^) and a density matrix ps,o G .^i(J^), i.e. such 
that Trps^o = 1 and ps,o > 0. Then we put 

O :=Os®W(V'>.), Po :=Ps,o®W(V'x)PaW*(0x). (1.13) 

If Ps.o = IV's,o) (V's,o| for some V's^o € J^s, then po = I^I^o) (*o| with vJ/q = Vs^o ^ >V(V'x )^ and this is assumed for 
simplicity in the next section (for notational convenience, we use the general case in later sections). We need to 
assume some regularity properties on V'x ^ 

Assumption C (Regularity of initial states and observables). 

dt{l + \t\r\h^,{t)\<^, zvith /^^(t) (0,e-'l9l*0^)^ (1.14) 

dt{l + \t\)"\h^{t)\<oo, with /i^(i) := ((^,e'''?'V>.)h (1-15) 



sup|/ix,(t)|(l +i)" < 00, with := (V^x,e'l«l*Vx)f) (1-16) 

t>o 



for a > 0. 



Throughout our paper, we always assume that Assumptions lAl and ICl are satisfied with the same parameter 
> 0. This is done for the sake of simplicity, though it slightly weakens the result (see Remark lLU belowl. 
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1.4 Results 

Write 4't := e^'*^^'o for some initial vector \I>o G ^ and H as defined in \1.5\ . We define the Weyl algebra 
QUq, a > to be the C* -algebra generated by 'atomic' observables A®1 and Weyl-operators 1 yV(V'x) with 
^» e f) satisfying dl.lSb . 

Theorem 1.1. Assume that Assumption lAl and Assumption\B\are satisfied. Then, there is a Xq > such that, for any 
coupling strength < | A| < Aq, the following holds true: 

1) There is a hounded linear functional O H> (O)oo on SIB a such that 

lim(*t,0*t) = (0)oo (1.17) 

i— j-oo 

for any initial vector e ^ with = 1 and O £ 2Uq. 

2) Let vjf^g = (K) 17 , 1 1 1 1 = I be the normalized ground state of the uncoupled Hamiltonian Hg + Hp and let O be 
of the form Ill.l3l l with lip^ satisfying jl.lSi , then 

3) IfO and (t/zflt /s; f/ze Vt< , f^^^ determine them) satisfy the three bounds of Assumption^ then 

|(0)oo-(*t,0*t)| <0(i-"), (1.19) 

Remark 1.2. One cowZd be tempted to interpret the functional {■)ao as the expectation in the ground state of the coupled 
system, but this is not correct since, for a < I, the coupled system does in general not have a ground state in (although 
a ground state does exist in the Hilbert space corresponding to a different representation of the operator algebra). On the 
other hand, if the system does admit a ground state \I'gs € M', then by choosing \l/o = \I'gs, l|1.17|l immediately implies that 

(O)oo = (*gs,0*gs). 

Our second result bounds the number of emitted bosons. Let N be number operator on the Fock space 

Theorem 1.3. Assume that Assumption \A\for some a > and Assumption^are satisfied. Then, there are Xq,Hq > 
such that, for any coupling strength A with < | A| < Aq, complex number k with \k\ < k'q, and initial vector ^Pq with ijj^ 
satisfying the bound (|1.14l l in Assumption^ we have 

|(*t,(l(8)e''^)*t)| < C'expl^Clip-™"^"'!"), t>0 (1.20) 

where the constant C depends on '^q, but C does not (and none of them depends on A, n or t). In particular, if a > 1, then 
the LHS is bounded uniformly in time. 

Remark 1.4. As indicated below Assumption^ we prefer to keep one constant a throughout the paper. Let us describe the 
possible improvement of Theorem 12.31 if we were to drop this constraint. Assume again that Assumption\A\holds for some 
a > and assume the the bound l|1.14b in Assumption^is satisfied for some > Q, then the photon number bound 
(|1.201 still holds with a determined by Assumption^ regardless of the value ofa^. This is clear from the inspection of the 
last part of the proof of Theorem [73] in Section \4.2.3\ 

1.5 Discussion of the results 
1.5.1 Quadratic Hamiltonians 

The easiest way to understand our results and the different conditions involved, is to compare them to an 
integrable model where the same questions can be asked. Consider the formal Hamiltonian 

f dq\q\a*ag + f dq (^^iq)a* + , e f) (1.21) 
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which fits into our framework by taking — C (in that case i/s is an irrelevant number). By completing the 
square (which, in this context, is a special case of a 'Bogoliubov transformation') we can rewrite it as 

H = Egs + l dq\q\blb„ b, := a, + Eg, ~ J dg^j^ (1-22) 

If (/> e (] and (j)/\/\q\ G f), then the term linear in a/a* is an infinitesimal perturbation of the quadratic term (the 
first term in (|1.21b ). The operator H is self adjoint on the domain of the quadratic term and it is boimded below 
by Egs > — oo. If moreover 4>/\q\ E f), then H has a normalizable ground state, given by 

vl,^, :=.e-ll'^/l«lll'w(<^/|<z|)r!, (1.23) 

with the Weyl operator yV(V^), ^p E t) as defined in Section [Ol We refer to e.g. [SI for an extended and rigorous 
discussion of quadratic Hamiltonians. 

Let us look into ergodic properties of the evolution generated by a quadratic Hamiltonian. Let V'x G L^, and 
consider the observable O ~ yV(V'x ). Then by explicit calculation 

(r!,e"^Oe-»^r!) =e-ll'^>'ll%2'^'=<'^-('="'"-^)w) (1.24) 

Recall the correlation function /i>^ (t) = (■0xi , e'*l^l(/)) from Assumption ICl Clearly, if h-^ e L^(]R, dt), then the RHS 
of J1.24|l converges as t oo. In particular, this can be true even when ct>/\q\ <^ \), that is, if H has no ground state. 
One can easily convince oneself that the t — > cx}-asymptotics of l ll.24|l does not change if we consider a general 
initial state po = l^oX^'ol with = W{i!K)^, ||*o|| = 1 and i/'x such that h^^t) = (0, e-^l'^l-f/'K) is integrable 
and htA{t) = (^/^x , e'^l'^'i/'x) vanishes at infinity. Next, we study the number of emitted photons 

(17,e"^e'^^e-'*^r!) = (^k j Ak{au + ^Ak))* {^k + ^t{k))^n), (1 - e'*l«l)|^ (1.25) 



exp(||<^t||2(e''-l)) =exp(^(e 



ds / ds'h{s - s') \ (1.26) 



It is clear that this expression remains bounded if (and only if) 4>/\q\ G L?, hence if H has a ground state. 
Moreover, we see that the rate of growth of the LHS of l ll.20ll in Theorem 1 1 .31 corresponds roughly to estimating 
h{s — s') by \h{s — s') \ in the integral. 



1.6 Plan of the proof 

Our two results, relaxation to the ground state and the photon bound, are very similar from the technical point 
of view, even if their physical meaning possibly is not. For this reason, we focus exclusively on the relaxation to 
the ground state in the present section, and we devote a few words to the photon bound at the end. 

The proof relies on the following philosophy. The original problem is formulated as a perturbation (with 
small parameter A) of an integrable Hamiltonian Ha, + H-p whose dynamics does not have the phenomenon that 
we want to exhibit: it does not relax into the ground state. This can already be seen by remarking that the atom 
S is not coupled to the field F and, as the former is finite-dimensional, its d3mamics is oscillatory. However, 
the Fermi Golden Rule dLlljl provides us with a picture that does capture the dissipative behavior: The fact 
that the state of the atom changes by jumps between eigenstates of the Hamiltonian suggests the following 
approximation. 

where M is the generator of a dissipative djmamics that we can loosely describe as the Markov jump process on 
eigenstates of with jump rates given by l|l.llb and exponential decay of the off-diagonal (in i/s-basis) part of 
the density matrix. The approximation becomes exact (at least as far as the S-state is concerned), as A -> 0, i — oo 
such that is held fixed. This was already advocated by Van Hove |19] and it was made precise by Da vies [6J. 
We state it explicitly in Proposition 13.31 and we review the proof in Appendix[Bl 

The underlying physical reason why l ll.27|l is a good approximation is that Pit is invariant under the free 
F-dynamics and any disturbance in the field caused by S is carried away (disperses) to spatial infinity quickly. 
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such that it is irrelevant for the further evolution, and one can pretend that the state of the field remains Psi • The 
dispersive property is a consequence of the temporal decay of field correlations for the uncoupled dynamics, 
which is our Assumption lAl It is therefore plausible that the evolution of S is Markovian on time scales longer 
than the time necessary for a field excitation to disperse away 

The approximate dynamics in l|1.27b exhibits relaxation to the projection onto the (uncoupled) ground state 
ipeo ^ ^ provided that sufficiently many jump rates are nonzero; this is captured by the Fermi Golden Rule 
Assumption IHl 

In a nutshell, our strategy is to use the dynamics (|1.27l l as a zero-order term of our expansion for the full 
dynamics. Note that our expansion is not simply in powers of the coupling constant A; the exponent in (|1.27l l 
clearly has zeroth and second order contributions in A. The reason we refer to (|1.27t as zeroth order is that all 
other contributions to the dynamics are small compared to this term, or rather, to its dissipative effect. 

Because of the jumps described by M, the dynamics (|1.27b is stochastic and hence our task reduces to con- 
trolling a small perturbation (the real dynamics at finite but small A) of a stochastic evolution. This is quite a 
tractable problem that can be handled by analytic perturbation theory of isolated eigenvalues and a cluster ex- 
pansion. Similar expansions were developed e.g. in I24ll20l and, very closely to the setup of the present paper, 
in [26,7]. 

The main result of this expansion is that we manage to represent Tr pt O as a one-dimensional polymer gas 
(this dimension corresponds to time). Then the problem of showing that there is a well-defined and unique 
asymptotic state is analogous to the problem of proving decay of correlations in the one-dimensional gas. 

More precisely, the polymer representation is 

Tr(ptO) = fcxfcx 51 n "(^)' = "/'^^ (1-28) 

A AeA 

where the sum is over collections ^ of sets Ac {0,l,...,n, n+1} such that for any two sets A-y.A^, dist(Ai, A2) > 
1 where we write dist(yli, A2) = minigAijeAa 1* — j|- The sets A are called poljnners, the numbers w( A) are poly- 
mer weights. Moreover, only the vi^A) with G A, (n + 1) e ^ depend on the initial state, resp. the observable. 
To a good approximation (the error made is not important for the present discussion), 

~ Trpo, fcx - Tr[(Peo ® Pn)0] (1.29) 

Our goal is to prove that 

lim Tr(ptO) =Tr(po)(0)oo (1.30) 

where (O)oo does not depend on po. Setting all v{A) = essentially amounts to pretending that pt = Pe„ <E) Pn 
for all t. In that case Tr(ptO) — kt^ky, and the expectation value of the observable is independent of the initial 
state, apart from the trivial normalization factor k^ . The polymers A contain corrections to this picture. These 
corrections originate from the fluctuations of the Markovian dynamics generated by M and from the corrections 
to the Markovian behavior I I1.27II . Pictorially, let A = {ri, T2, . . . t„j} with < t^+i, then v{A) describes correl- 
ated deviations from the assumption that pt — Peg ^ Pn in the time-interval Dom(ri), . . . , Doni(T,„_i), where 
Dom(T„i) = (1/A^)[r„i — 1, t,„]. The detailed construction of the polymer representation l|1.28|l is carried through 
in Section|2l 

To prove decay of correlations, we use a standard cluster expansion, which we review in Appendix lAl 
A possible condition for the applicability of the cluster expansion method is the "Kotecky-Preiss" criterion. 
Applied to our model, it demands that, 

«(A)e'^(^) < a{A'), (1.31) 

A:dist{A.A')<l 

with A, A' C {1, . . . ,n} and a(-) an n-independent positive function on polymers. We formulate (|1.31t (together 
with some other statements) in Lemma [2.41 In our case a{A) = eC\A\ where e = |^|2min(Q,i) seen as a 

renormalized coupling constant. 

To prove that the Kotecky-Preiss criterion l ll.31t is satisfied, we use a Dyson (or Duhamel) expansion, it relies 
on the smallness of the coupling constant A and the decay of field correlations. This is intuitive; if the correlations 
of the free field decay roughly as (cfr. Assumption [A) then one could conjecture, for example, that for 

A = {ti,T2} also the weight v{A) decays as (dist(Dom(Ti). Dom(r2)))~^^+"-' as T2 — ti — >■ 00. This picture 
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turns out to be essentially correct. In fact, we get |f ({ti, t2})| ^ Ce(T2 — ti)"^^"*""-*. This analysis (proof of the 
Kotecky-Preiss criterion) is accomplished in Section |3A] 

In the concluding Section |4l we prove our results. The main point of this section is to pinpoint how the 
cluster expansion gives rise to decay of correlations. Such reasoning is completely standard in high-temperature 
expansions of statistical physics, see e.g. |27J. 

To prove the photon bound, we use an analogous approach but this time we develop a polymer representa- 
tion like (11.281 1 for the quantity Tr(pte'*^). In fact, in our proofs, we provide one general polymer representation 
for Tr(Oe(''/2)^pte(''/2)JV) ^nd then we set k = to study Tr(ptO) and O = 1 to study Tiiptc'^^). 

In the rest of the paper, we assume that Assumptions lAl iBl and ICl are satisfied for a > and all the Lemmata 
will depend on this parameter a. 
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2 Polymer Representation 



In this section, we complete the first important step of our proof, namely we rewrite all quantities of interest 
through a polymer representation. First, let us discretize time by introducing a mesoscopic time scale A^^, where 
A > is the coupling strength. That is, we consider times of the form t = nj}? with ri e N (the discretization 
will be easily removed at the end of the argument). The main quantity that we study is 



Zn = ^„(0, poi't) := Tr 



Qg(K/2)Arg-i(n/A^)ff^^gi(n/A2)ffg(K/2)W 



(2.1) 



where N is the number operator on the Fock space and k e C is a sufficiently small parameter. The operators 
Po and O are the initial states, respectively observable constructed as in Section [L3l 

As announced, we develop a polymer representation for Z„. In fact, we will first construct a polymer 
representation with operator-valued polymer weights. In this representation, the polymers correspond to deviations 
from Markovian behavior. This is described in Section lZT] Bounds on the operator-valued polymer weights are 
stated in l2.2l but their proof is deferred to Section|3] 

Then we define the 'true' (scalar) polymer representation for Z„ by adding the excitations of the Markovian 
approximation to the already defined operator-valued polymers. This is done in Section l23l Why this leads to 
scalar polymers will be explained there. Finally, we need to provide bounds on the scalar polymers in order 
to satisfy the Kotecky-Preiss criterion. These bounds follow from the bounds on the operator valued poly- 
mer weights and this is described in Section |Z41 In all of the following sections, we prefer to treat 'bulk' and 
'boundary'- polymers separately, where the terms 'bulk' and 'boundary' refer to the time-dimension. The ana- 
lysis of the bulk polymers is the more subtle piece of work but the treatment of the boundary polymers demands 
additional notation. Therefore we treat the former first, and then indicate the (in all cases quite minimal) changes 
necessary for the latter. 



2.1 Operator valued polymer model 

2.1.1 Definition of the deformed dynamics 

Starting from ll2.1b , we would like to move the deformation parameter k into the dynamics e~'*^. Given an 
operator X, we introduce the left, right and two-sided multiplication, acting on operators p 

C{X)p:^Xp, n{X)p:=pX, M{X)p:= XpX. (2.2) 
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Defining the Liouvillian L = ad(iJ) = C{H) — TZ{H) with H as in l|1.51 , we may rewrite jl.l^ as 



Z„ =Tr OX(e('"/2)iV)g-i(n/A )L^^ _ (2.3) 

We formally define the deformed (not self-adjoint) Hamiltonian 

- e('=/2)^i?e-(''/2)^ (2.4) 

and the deformed Liouvillian 

:=£(ii',)-7e(if_,) =X(e(«/2)^)i:X(e-(''/2)^). (2.5) 

Then, formally 

7\/[(e(«/2)W)e-i*iM(e-(-/2)^) =e-'*^» (2.6) 

which would allow to rewrite 



Zr, = Tr 



(2.7) 



To make these manipulations precise note first that relying on the relative boundedness of e'^/^^^ffie^^''/^^^ 
w.r.t. Hs + Hp, which follows from 1II.6I I, it is easy to construct H^^ and as unbounded closed operators and 
to show that they form an analytic family (in k) of class A in the sense of Kato |21J. Then (|2.7b is a simple 
consequence of the functional calculus in case k £ iR. The upcoming Lemma [2.11 gives a constructive meaning 
to it for arbitrary k and provides an expansion that we will use in practice. To state this lemma, we need a few 
additional definitions. First, some additional Liouvillians: 



and 

Defining 



Ls := ad{Hs), Lp ~ adiRp), Lj := Xad{Hi) (2.8) 
(e'^l9l+''/V(g)a; + e-'^l«l-«/2VK^aJ (2.10) 



we have 



Li,,(s) - \C (e-^^ (D) ® $,(0, s)) - n (e-^^(i?) ® s)) . (2.11) 



For * € we decompose * = EmeN*™ where = l[N=m]'^, i-e. € ^ » with J^F^m C ^ 

the m-boson sector, see e.g. [9] for more details. Define the dense subspace &i{J^) c to be the space 

consisting of finite linear combinations of rank-one operators satisfying 

Vto e N : max(||^'„||,||*;„||) < for some C > (2.12) 

V m! 



Lemma 2.1. The LHS of equation (|2.6b defines an operator 

Q-itL^ . ^^(^) ^ ^i(^) 

/or all K £ C and t e M. Groen p g ^i(J^) t/ie map k 6"'*^"^ fs holomorphic from C to ^i(J^). Moreover on this 
domain 

g-i(ti+t2)L„ ^g-itiL«g-it2L. (2.13) 

e-"^V = e-'*(^^+^-) ^(-i)™ y dii...di„,Li,,(t™)...Li,Jt2)ii,.(ii)p (2.14) 

0<ti<. ..<t„<t 

w/zere f/ze to = 0-term on the RHS is understood to equal p and sums and integrals converge absolutely. Finally, the RHS 
ofeq. \1.7\ is well defined and Il2.7b holds. 
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Proof. Starting from = Lg + Lp + XLi,^, and iterating the Duhamel formula 

e'»^''p = e-"(^''+-^^V-i / dse~'(*~")^'^Li,«,e"'"(^^+^=^V (2-15) 

we formally arrive at l|2.14b . Hence, the only nontrivial claim in the lemma is the absolute convergence of the 
series on the RHS of l|2.14b and the fact that it belongs to We refer to |10J for an explicit proof, which 

relies exclusively on the well-known estimate 

j Aq^^{q)a*-^,^ < V^^^ ||*m|L^ ^ dglV^g)!', fora# = a,a* (2.16) 



□ 



Hence, in what follows, we freely use the operators e and the group property II2.13II . 



2.1.2 Splitting of the dynamics 

We define the reduced dynamics Qt : S§i (Mi) (Mi) of the atom; 

Qtps ■■= TrF[e-"^''(ps Pn)] (2.17) 

where Tip : Bi (J^/f ) — > Bi {J^/fs ) is the partial trace and the well-definedness of the RHS follows from Lemma IZTl 
A large part of our analysis serves to prove that Qt tends to a one-dimensional projection as t ^ cx). 
We start by rewriting Z„ in l|2.7b . Recall Lp = ad(iJF) and introduce operators Ur with r e N. 

JJ^ ^ gi(r/A2)LFg-i(l/A^)L„g-i((T-l)/A2)LF (2.18) 

The motivation for this definition is that the product of Ur telescopes into 

[/„... L/i = e'^n/\^)L^^-i{n/\')L^ (2.19) 
In particular, if we choose po = ps,o ^ Pn and O = Os "8) 1 then eq. i2.7^ becomes 

= Trs[OsQ«/A2ps.o] (2.20) 

and hence in this case the study of Z„ reduces to the study of Qt- The main idea of our approach is that, at least 
qualitatively, the main contribution to Qn/x'^ can be inferred by approximating Qn/\2 by (Qi/a^)"- We rename 
T := Qi/\2 and we define the 'excitation operators' 

Br = Ur - T ®1 (2.21) 

Our task is to understand how the behavior of T" is modified by the excitation operators Br ■ To quantify the 
influence of the latter, we now develop a formalism. 

2.1.3 Correlation functions of excitations 

We abbreviate 

= ^(^i(^s)), = ^(^i(^f)) (2.22) 

Define, for W, W e^s'^ the object 

W (E)s W e <E) <S) 

as an operator product in F-part and tensor product in S-part. Concretely, let W — Ws ® Wp and W = VFg Wp. 
Then 

W (g>s W := Ws(E>Ws(g) WfWp. 
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and we extend this by linearity to arbitrary W, W'. Iterating this construction we define for Wt E Ms (g) My, 
i = 1, ... ,171 

(Note that no analysis problems arise since is finite-dimensional.) We define the 'expectation' 

E : (^s)®" ® S^v ^ (^s)®" 

as 

E(VK)J :-TrF[W^(J®Pn)], Je (^i(jrs))®" 

Obviously, the action of E is extended to unbounded W satisfying ® Pq) € ® J^f)- 

In particular, by (an obvious generalization of) Lemma l2.ll we can consider W = Ur^ (X)s . . . <8)s for any 
TO-tuple of times ri, . . . , r^. For example, we rewrite the definition of T as T — E([/^) (this is with m — 1). Let 
A = {ri, r2, . . . , Tm] C N with the convention that t; < r^+i and define the 'time-ordered correlation function' 

Ga := E {Br^ ®s B,„_, (»s • • • ®s BrJ e (^s)®" (2.23) 

Note that Ga = when the set A is a singleton, as follows directly from Br = Ur - E(C/t-) ^ E(i?T-) = 0, and 
Ga = Ga+t because e-'*^^Pn = Pn- 

It will be convenient to label the ^s's and to drop the subscript S (since we will rarely need ^f)- Therefore, 
let T e N be copies of and let A c N be as above. We define by 

3£a ■■= ^T„_i «)■••«) ■ 

Obviously, is naturally isomorphic to by identifying the right-most tensor factor to ^i, the next one 
to etc. We denote this isomorphism from ^® to Sf^A by I a and we will from now on write Ga to denote 
/a[Ga] since Ga acting on the unlabeled space ^® will not be used. Similarly we will most often abbreviate 
/^[V^] {y £ by Vt, which will lead to a slight abuse of notation, see below and in Section [2.1.61 Consider 
a collection A of disjoint sets A , then each of the spaces ^AeA is naturally embedded into Msu-ppA/ where 
Suppyl — Ua^a^- Consequently, given a collection of operators Ka & ^a, A e A, we can define 

(S)AeA Ka e ^Supp.A (2.24) 

In particular, we have Ma = ®T(^AMr. In the literature on quantum lattice systems, where similar constructions 
are necessary, one usually identifies Mr with the subspace ...®\®\®M®\®\®...oi the infinite tensor 
product, such that, for example, ®AeAK-A is simply written as Ka- However, we chose to have the tensor 
products explicit in the notation. 

We define the "contraction operator" T : Ma M, by first giving its action on elementary tensors. Consider 
a family of operators Vr e M, and set 



r 



® Vr 
reA 



^r,„ Vr^_, ...Vr,, where > t„_i > . . . > n (2.25) 



and then extend linearly to the whole of Ma- On the LHS, we abbreviated Ir [Vr] by Vr- 

By expanding Ur = T ® 1 + Br for every t e /„ := {1, . . . , n} in the expression for the reduced dynamics 
(|2.17t , we arrive at 



Ga ® Tr 



(2.26) 



Note that, for each operator appearing in the tensor product, we have specified the space Ma on which it acts. 
In contrast, the order in which we write the factors inside the T[ ] does not have any meaning. 

Remark 2.2. The above construction with tensor product spaces Ma and the contraction operator T does of course not 
depend on the fact that the spaces are indexed by elements ofN. The only requirement is that the index set is ordered. In 
particular, in Section \3?l\ we will use the same formalism, but now with copies of M indexed by (a finite number of) times 
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2.1.4 Polymer expansion for cumulants 



The cumulants or "connected correlation functions", denoted by GJ^, are defined to be operators in !Ma satisfying 
'^{r} = G{r} = and 

Ga'^ J2 f ® G'X] (2.27) 

partitions A of A 

The tensor product in this formula makes sense since ^a' = ®a^aS^a whenever is a partition of A! (cfr. 
previous section). Note that this definition of connected correlation functions reduces to the usual probabilistic 
definition when all operators that appear are numbers and the tensor product can be replaced by multiplication. 
Just as in the probabilistic case, the relations l|2.27|l for all sets A fix the operators G\ imiquely since the formula 
(|2.27l l can be inverted. 

Plugging l|Z27t into llZ26l l, we get 

where *8° is the set of disjoint collections of subsets of /„. The formula Il2.28b is the starting point of our analysis. 
Note that in the case where A contains at least one set A that is not a discrete interval (a consecutive set of in- 
tegers), there is no obvious way to write the RHS of Il2.28b as an operator product. This was the main motivation 
for introducing the formalism with tensor products and the contraction T ■ Following a standard terminology 
in statistical mechanics we call the sets A c /„ polymers and the function weight of the polymer. Unlike in 
statistical mechanics this weight is operator-valued and our objective is to manipulate this expansion to arrive 
to standard scalar valued weights. 



Af^A 



G", 



rG/„\Supp^ 



(2.28) 



2.1.5 Notation for combinatorics 

We gather some notation (partially already introduced above) that will be used throughout. We let /„ = 

{1,. . . ,n} which is interpreted as the set of discrete times, often denoted by T,r' It will be convenient 

to add two "boundary elements" to this set, which are represented by the times and n + 1, hence we set 

/„ = /„ U {0, n + 1} (2.29) 

The set ?8„ is the set of collections of subsets of /„, i.e. *8„ = 2^^'"\ Elements of *B„ will mostly be denoted 
by ^ = {Ai, ^2, . . . , Am} with Ai c /„. We define two important subsets of *8„: *B° is the set of collections of 
mutually disjoint sets in /„, and *B^j is the set of collections of sets such that the distance between any of them is 
greater than 1; 

"Bi ~ {A C 2^" \ yA,A' eA: {A^ A' ^ dist(A, A') > j)}, j = 0, 1 (2.30) 

where dist(j4, A') = mmr£A,T'GA' \t — t'\. Similarly, we define *B„, *8" , starting from /„ instead of /„. Inter- 
vals in /„ and /„ are sets of consecutive numbers. For any A, we say that the intervals Ji . . . , Jm are the maximal 
intervals in A iff. UjJj = A and the collection { Ji, . . . , J,,,} is in 03 ,\ or *8;',. 
We define in general the support of ^ € *B„ as 

SuppA = UaeaA (2.31) 

and the 'span' of sets as 

SpanA :— {mm A, . . . , max A}, Span^ — Span(Suppyl) (2.32) 

that is SpanA, Span^l is the smallest interval that contains A, Supp^l, respectively. The size of SpanA, Suppyl is 
called the diameter of A, A, denoted by 

d{A) 1 + max A — min A, d{A) = 1 + maxSupp.4 — minSupp.4 (2.33) 
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2.1.6 Initial state and observable 

To deal with the observable O and initial state povna convenient way we also define the operators C/q, Un+i 

Uop := e(«/2)^W(Vx)pW*(V'K)e("/')^, (2.34) 
Un+iP W(e'("/^')l«l^>,)p (2.35) 

where we wrote W{iIj) instead of 1 ^{ip). Note that Un+i depends on the total macroscopic time n, which is 
a notational drawback of our formalism. We introduce a modified reduced dynamics Q„, as 

Q„ps Trp [c/„+ie'("/^')^-e-'("/^')^''f/o(ps Pn)] (2.36) 

such that we have 

Zn{0,pQ,K) = Trs[OsQriPs,o] (2.37) 

which generalizes l|2.20b and reduces to the latter when = V'xi =0. To check l|2.37b , note that H'(e'*I^IVxi) = 
e"^''W(V'xi)e"'*^'='. 

It is straightforward to extend the formalism of Sections I2.1.3H2.1.4I to incorporate the times 0,n + 1. We 
define 

To E(t/o) = (e(''/2)Ar^(^^)^^g(«/2)A'-^(^j^^j (2.38) 

T„+i EiUn+i) = {n,W{ijy,)n) 1 (2.39) 

and set 

Br := Ur-Tr, T = 0,n + 1. (2.40) 

Next, we define copies ^o,^n+i of and also for A n {0,n + 1} ^ ^, completely analogous to the 
construction in Section [2.1.31 such that the definition of (connected) correlation functions extends to all 
A c /„. In our expansion, we will need 7o[7o], -fn+iiTn+i] which we write simply as ro,T„+i. With these 
definitions, the representation (12.281 1 is generalized to 



AeA Te/„\Supp^ 



(2.41) 



2.1.7 Norms 

Let us introduce some conventions. For S acting on J^, we write 



IS*!! sup II^V-II (2.42) 

V-e-^sJlV-ll^i 



and for E acting on ), we write 



\\E\\:= sup \\E{p)\\i, with IIpIIi = Tr|p|, (2.43) 

pe.®i(.^s),||p||i=i 

i.e. the natural operator norm on ^(^i(J^)). 

For E e \ A\ > 1, we exploit that E can be written as a finite sum of elementary tensors , i.e. 

S = E^^ E)reAE^,r, E^^r G ^r, (2.44) 

1/ 

to define 

TEA 

where the infimum ranges over all such elementary tensor-representations of E. This norm is useful because of 
the following properties (trivial from the definition): 
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1) For any family of operators Kaea with Ka € and A a collection of disjoint sets, i.e. A E S", we have 



«) Ka 
AeA 



< n ll^^ll- (2.46) 
AeA 



2) For any Ka e ^a, 

\\T[KA]\l<\\KAh (2.47) 

2.2 Bounds on operator-valued polymers 

We aim to set up a perturbative scheme where the will describe small corrections to T, and hence we must 
specify in what sense the operators are small. 

We will often need to distinguish between bulk polymers, i.e. subsets of /„, and boundary polymers, i.e. 
subsets A C In for which A n {0, n + 1} 7^ 0. Lemma [2.31 gathers the necessary properties of operator valued 
polymers. To relate these properties to assumptions on the original model, we introduce the 'renormalized 
coupling constants' 

g _ |^|2min(a4)^ £ ;= max( | A| , e) (2.48) 

We use in general G, c to denote constants that can depend on all model parameters except the coupling constant 
A, conjugation parameter k, the macroscopic time n and the initial state and observable. By G, c we denote 
constants that can also depend on the initial state and the observable (but not on A, k or the macroscopic time n). 
For operators J, J' e ^(Jfs), we use the Hilbert-Schmidt scalar product (J, J') = Tr[J* J']. Recall the diameter 

of A, d{A), and write d{A)" := (d(A))". 

Lemma 2.3. For sufficiently small |A|, but \ ^ 0, the following hold uniformly in A, k (i.e. gT and all constants can 
be chosen independent of \,k) 

1) The operator T has a simple eigenvalue e^^ with 9t — Ot{^, cind a gap gT > 0, in the sense that 

||r™ - e™^^i?|| < Gre^d^^l-s^), to e N (2.49) 

where R = R{\, n) is the (one-dimensional) spectral projector associated to the eigenvalue e^^. Moreover, Ori^, n = 
0) = and liniA^o 6't(A, k) = 0. R is of the form 

R=\V){V\- (2.50) 

where, for k = 0,we can choose fj = 1 and rj a density matrix. For k e R, zve can choose ry, r) positive. 

2) Let e, e be as in Il2.48b . The bulk polymers satisfy 

max (C^f)-^'^'~'^rf(^)"llG^IU < 1 (2.51) 

" Adln-A^T 

For boundary polymers, take r = 0, n + 1, then, 

J2 (G6)-l-^"^"ld(A)"||G^|U<G (2.52) 

Aci„:A3T,A^{a,n+l} 

For A = {0, n + 1} (excluded from the sum above), we have d(A)"||G^||o < G. 
The proof of this lemma is in Section |3l 

2.3 Scalar polymers 

The representation (|2.28b evokes the picture of a leading dynamics T interrupted by excitations, indexed by the 
sets A G A. One could call this representation a polymer expansion, but it is not yet what we want because the 
values of the polymers, i.e. G^, are operators. To make them scalar, we exploit the dissipativity of the model. 
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namely the fact (see Lemma [2.3b that the reduced dynamics T has a maximal (in modulus), simple eigenvalue. 
Recalling that the corresponding spectral projection is denoted by R we write 



T = RT + Rj_T, RT^e^'^R 
where Rj_ := 1 ~ R and we have RT = TR. Hence, in particular 

RTRaT = R±_TRT = 0. 



(2.53) 



(2.54) 



We insert this decomposition into the expansion | |2.28|I : Let J {A) ~ { Ji, . . . , Jfe} be the family of maximal 
intervals in /„ \ Supp.4. The relation Il2.54b implies that each tensor product ® Tr in 112.281 1 may be replaced by 

® {RT)r+ ® (R. I T)^. Thus fTM becomes 



Ae<30 JGJ{A) 



(E) Ca (E) {R±T)r (E) {RT), 

AeA reSuppj Te/„\(SuppjuSupp^) 



(2.55) 



We will now produce a new family of polymers A' = A {A, J) by "fusing" some of the sets in the family AU J'. 
Let r = T{A, J) be the graph with vertex set V(r) = AU J and edges {S, S'} whenever the sets S and S' are 
adjacent i.e. dist(S', S") = 1 (in particular this implies that at least one of them is in A). For each connected 
component 7 of F write V(7) = U J-y, set A'^^ SnppA-y U SuppJ'-y and let A be the family of A'^. Defining 



ViiA^,Jy)) 



E> E) {R±T)r, 

A&Aj reSuppJ"^ 



the identity l|2.55b becomes 



Ae<s^ JcJ{A) 



E>V{{A.y,J-y)) E> {RT)r 

7 TGl„\SuppA'{A,J) 



(2.56) 



(2.57) 



The next step is to write II2.57I I as a sum of families A, i.e. to fix A and sum over A and J'. 

Definition 2.1. Let ^ be a family of mutually disjoint subsets of /„ and let J' he a family of mutually disjoint 
and non-adjacent intervals in /„ and disjoint from the sets in A, i.e. Suppy^ n Supp J' = 0. We say that S = {A, J) 
is a fusion iff. 

1) dist(/„ \ Supp5, Supp J') > 1 where SuppiS := Supp^ U SuppJ'. (Fictorially; the intervals J are in the 
"interior" of Supp5.) 

2) The graph T{A, J) is connected. 
The set of fusions is denoted by (3/^. 



Defining 

5=(AJ)6S{:Supp5=A' 

we obtain the representation for Qnjw^ in terms of fusions as (we drop the prime from A and A) 



(2.58) 



A<i_<&\ 



E> {RT)r E> 
Tei„\SuppA AeA 



(2.59) 



Note that by construction A £ 05 since the sets A'^ above are non-adjacent i.e. their mutual distances are at least 
2. Hence, for any A in the formula above, all r that are adjacent to Supp^l carry the projector RT. A pictorial 
way to phrase this is that any A in l |2.59l l is surroimded by projections R. We exploit this by defining 



v{A') 



T 



^ViA') (g) Rr 
Tei^\A' 



v{A') e ^ 



(2.60) 
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Note that v{A') is a multiple of R unless A' 3 1 and/or A' 3 n. To eliminate these boundary effects, we consider 
^Qn/A2^ instead of Qn/x^, for the time being. Then Il2.59b implies 

i?Q„/A2i?= ^ (i?r)("-|S"PP-^l) ]^ (2.61) 

where the product on the RHS is commutative since the projection R is one-dimensional. It is convenient to 
extract the contribution due to the RT, by defining (recall ?/. fj from l|2.50b ) 

= e-l^l^^(?7,{)(A)r/) (2.62) 

Indeed, since (?/, rj) = 1 (because RR = R), we get 

RQ^/x^-R - i?e"^- n (2.63) 

Ae'Bl AeA 

and the RHS is indeed a 'scalar polymer representation', i.e. the weights v{A) are numbers, and we can study it 
with the help of a cluster expansion. Note for later use that 

Tri?g„/;,2i?= Z„(r)«) l,?7(g)Po,K). (2.64) 

Our real object of interest however is the partition function ^237) expressed in terms of the operator (|2.4H . We 
proceed with the latter as with Qn/x^ ■ Decompose as in (12.531 1 

Tr = RTr + R±Tr, T G {0, n + 1} (2.65) 

where R is the same operator as in Il2.53t . Since ro,T„+i are proportional to the identity operators, we have 
again RTr = T^R. As a consequence, for all r, r' e /„ 

RTrRl_Tr' = Rl_TrRTr' = 0. (2.66) 

We end up with the expansion like (|2.59t for Qn with the expected difference that we have to replace /„ by /„ 
and by *B,\ (also in the definition of fusions). 

For Zn of (12.3711 we get the scalar valued polymer expansion as in (I2.63t with the change that the weights 
v{A) for An {0,n + 1} ^ have to be slightly modified because the weight of the polymer is influenced by the 
initial state and the observable. We keep the definition of v{A) given in | |2.60|I . Then, recall that R = \ri) {fj] and 
define 

{{fi,v{A)ps.o) OeA,n + l(^A 

{Os,v{A)f]) n+leA,0(^A (2.67) 

{Os,v{A)ps,o) {0,n+l}eA 



where 
and 



kA := (l[OeA]fcK + 1[0^A]) (l[n+lGAl^x + l[n+l^Al) (2.68) 

(77,roPs.o) =Trs(77,Ps.o)(e("/'^^W(V><)f^,e(''/2)^W(^K)f^) (2.69) 
ky, := {Os,Tn+iv)=TTs{Osv)T^M^{MPn) (2.70) 

With these definitions, we arrive at 

which is our representation of Zn as a polymer model. A few remarks: 

1) The case where the factors ky ,ky, vanish, will be of no concern, as we will explain later, following Lemma 
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2) Up to now, we did not indicate the macroscopic time n in the notation. Let us, only in the next Hnes, 
indicate this dependence by a superscript. Consider A such that max A — 1 < n < n', then 

(2.72) 

Furthermore, v^^^A + t) = v^^^A) for t e N, as long as < minA, maxA + r — 1 < n. If we allow 
polymers to contain the final time n + 1 but not the initial time 0, i.e. < min A, max A — 1 < n < n', then 
we have 

= + (n' - n)). (2.73) 

These properties follow from the property Ga = Ga+t and the expression for the 'final-time' operators 

Un+l, Un' + l- 

3) If we choose po = i] (g) Pn, then v{A) — whenever E A. Indeed, Indeed, in that case we can substitute 1 
for Uq, hence Bq — 0, and we have R±ps,o — because ps,o = V- By analogous reasoning, we check that if 
we choose O = f] (E) 1, then v{A) — whenever n + 1 E A. Note that these two observations are consistent 
with (|2. 6312. 64b which tell us that one can omit the boundary polymers. 

We finish this section by a pictorial illustration of the construction. Each term in the representation of Qn/x^ 
in Il2.28b can be represented by a picture as in The horizontal axis is the time-axis which has been divided 
into intervals of the form A^^ [r — 1, r] by the discretization procedure. Each term in Il2.28b is specified choosing 
a collection A, this is indicated on our picture by connecting the r that belong to the same A E Ahy straight 

diagonal lines above the axis. The times r E Supp^ are drawn with the symbol ■ and the r ^ Supp^ are 

drawn with a 

As to the operator value of this picture, each -^^^ corresponds to T, and each coimected component of • 

corresponds to an operator G^. 




We now generate new pictures by splitting 

T = TR + TRj_ = = W/W + (2.74) 

In this way we obtain 2^ pictures from Figure [TJ with £ the number of One of these is the upper picture 

in Figured Note that splittings in which one -W/m is adjacent to an * do not contribute because RR± = 0. 
Therefore, we are in fact choosing the splitting for each maximal interval in /„ \ Supp^. In a next step, we do 
not distinguish between excitations that originate from Br or TR± and we write simply 

» for both t and • . Moreover, we fuse adjacent ♦ to form new polymers A', the lower picture 

in Figure|2]shows the result of fusing the upper picture. As a result, these new polymers are surrounded by waw, 
and since those correspond to one-dimensional projectors, the operator-valued contribution of a new polymer 
to the total sum is independent of the presence of any other polymers. This is why the new representation is 
scalar. 
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Figure 2: In the upper picture, we split T = TR + TR±^ for all t ^ Ul^^Ai and we choose one of the two terms 
for each r. In the example displayed, we choose TR±^ for ti_2,3 = 14, 15, 21 (those that correspond to the triple 
lines in the upper picture), and TR for all other r's. We fuse blocks marked with TR±^ with a polymer whenever 
they are adjacent to that polymer and we fuse two polymers whenever they are adjacent (possibly after having 
been fused with Ti?_L -blocks). In the example above, A'^ = Ai U A2, = A3 U A4 U {n, r2, ra}. 



2.4 Bounds on scalar polymers 

To obtain bounds on the new, scalar, polymers, we use the smallness of the operators G% expressed by the 
renormalized couplings e, e, and their summability in A when one of the elements of A is held fixed. The TR±^ 
factors, which glue together several GJ^ -operators into one scalar polymer, are not small (they are of order 1) but 
whenever we line up i of them, their weight decays exponentially in £ by the ergodicity of T (see Statement 1) of 
Lemma |23). Save for 2 special cases, dealt with in \2.77\ below, all scalar polymers contain at least one operator 
G\ so that one can always extract a factor e or e. 

Lemma 2.4 (bounds on polymer weights). For sujficiently small \ \\ 7^ 0, there is a ay > such that 

1) For bulk polymers, 

J2 e'^"l^ld(A)"|w(A)| < Ce, (2.75) 

AC/„:A3r 

2) For boundary polymers, take r = 0, n + 1, then 

e''"l^ld(A)"|w(A)| < Ce, (2.76) 

A C In ■■ A 5 T 

Aj^{0,n + l},Aj^I„ 
and the boundary polymers that are excluded in the sum above satisfy 

Hin)\ < Ce-^"^^", \v{{0,n + l})\ < Cn-^ {1.77) 
Proof. We first prove (12.751 1. Choose A' e /„. By the properties of the norm || • ||o stated in Section l2.1.7l , we get 
HA')\ < Ce-I^l^°^- J2 Pll^nil^AlUniK™^)'^'!!- (2.78) 

S={A,J)eei:SuppS=A' A£A J&J 
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where p is the number of maximal intervals in SpanA' \ A', and we will upper bound p < C|Suppy^|. For any 
fusion S — {A, J) such that Supp(^ ij J) — A', we have 



d(Ar < n '^(^)" 11(1 + 1^1)" (2.79) 
AeA JeJ 

which follows immediately from the connectedness of the graph r(^, J') in Definition 12. II 

To continue, let us interpret pairs S — {A, J) as a collection of sets S ^ Ay^J together with a label indicating 
whether they are 'A-sets' (g Ai) or ' J-sets' (e (of course, only intervals can be J-sets). Let us gather such 
collections of labelled sets in the set 6„. Furthermore, we consider the adjacency relation ^ defined by S* ^ 5' 
iff. d(S, S') ~ 1 and at least one of the sets S*, S' is an A-set. Then, it immediately follows that, for any fusion S, 
the graph (S,^) is connected. 

We introduce (with , Cw > 0, to be fixed below) 

f d(A)"(eC^)-(l^l-i)||G".|L if is an A-set 
w,iS) ■■=< (2.80) 
[ C-ie-^-l-'l ifSisaJ-set 

Then we can bound the LHS of (I2.75t as 

e^^A'\d{A')°\v{A')\<eC ^ H ^^('^) (2.81) 

5Ge{:SuppS=yi' ses 

provided one chooses + < gr and e small enough (depending on and c^). To check (|2.81l l, use II2.781 
12.791 1, the bounds |l-R||o < C and \ J\ < C|Supp^| for any fusion {A, J), the fact that any fusion contains at least 
one A-set, and the bound from Lemma [2. 31 1); 

g-|j|ReeT||(j.^^)|./||| < c-re-l^le^. (2.82) 
The desired bound (12.75b will trivially follow once we establish 

1[(5,~) connected] Wt{S) < C (2.83) 
56e„:Supp59T SeS 

because (5, is connected for any fusion S. To prove (12.83b , we first establish that, provided is sufficiently 
large, 

, /Q\ ,^ , s \ 1^1 if 5 is an A-set 

J2w,{S)e'^^^^ <a{S'), a(5):^ ' (2.84) 

1 if 6 is a J-set 

Indeed, if 5 is an A-set, one uses l|2.51b and if it is a J-set, one sums the exponential. Relying on l|2.84| l (which is a 
'Kotecky-Preiss' criterion) in the terminology of Appendix|Al we now apply the general combinatorial estimate 
Lemma Fa. 11 obtaining 

J2 l[(SU{S'},'-) connected] H ^^(^) ^ ^"^^'^ (2.85) 

5Ge„ ses 
(we bounded the indicator 1[ j by k{S U {S'}) and we used (|A4b ). Consequently, 

5 G 6„ : Supp5 9 T ses S':S'3T 5 g 

(5, ~) connected (5 U {S'}, ~) connected 



< ^ u;e(5")e"(^') < C (2.86) 

S':S'9t 

This proves Il2.83b and hence also (12.75b . 

The other claims of Lemma|23]are proven by analogous reasoning, using the bound on boundary polymers 
(12.52b . Indeed, the general idea of the proof of (12.75b was to define in Il2.80b such that it satisfies Il2.83b and 
to extract a factor from all terms in the sum on the LHS of Il2.83b , which was possible because all fusions S 
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contain at least one ^-set. To get \2.76\ , we modify the definition of for S such that S'n{O,n+l}7^0to (note 
a new boundary-dependent constant Cw) 

( d(A)'^(C^A-^(eC^)-\^'^'"\\\G'\\\o if is an ^-set 
w,{S) I " (2.87) 

[ C-ie-'^-'l-'l if 5 is a J-set 

and we use S„, 6^ instead of &n,&l by replacing in the definitions /„ by /„. Then (12.81b still holds with e 
replaced by e provided that S contains at least one A-set with | A n J„ | > 0. This is the case unless A' = Supp5 is 
{0, n + 1} or /„, which are the two special cases in Lemma [2.41 Next, we need to establish Il2.83b with C on the 
RHS replaced by C and r = 0, n + 1. It is pedagogical to split this estimates according to which boundary times 
S contains. Let us consider the case where G Supp5, but n + 1 ^ Supp5, in particular choose t = 0. Then, we 
write, analogous to (|2.86b , 

5 G Sn : n + 1 ^ Supp5 ■5e<S S':S'BO 5 £ e„ SeS 

(5, ~) connected , e Supp5 (5 U {5'}, connected 

The sum over S e (5„ concerns only bulk quantities, and it is therefore bounded by e"^'"^ -', just as before. We 
then conclude by bounding 

E we(5")e''(^'^ < C (2.88) 

S':S'30 

as follows from (|2.52l l when S is an A-set, otherwise it is the same bound as before. The cases where S contains 
n + 1 but not 0, or n + 1 and 0, are dealt with analogously, and we get l|2.76b . Finally, we turn to (|2.77ll . The 
claim about A = {0, n + 1} follows directly from the bound on that polymer in Lemma IZ3l The claim concerning 
A = In differs from Il2.76b because it is now possible that the fusion S consists of a single J-set J = /„, in which 
case one cannot extract e. Taking this into account, this claim follows as above. □ 



3 From Hamiltonian dynamics to polymer models 

In Section|2l the reasoning was largely independent of the details of the underlying Hamiltonian model. Indeed, 
in that section, we assumed some properties of the operator-valued polymers in Lemma [2.31 and we explored 
the consequences of these properties. Now, time has come to prove Lemma [Z3l This lemma discusses the er- 
godicity properties of the operator T and bounds on the correlation functions G^. In both cases, the proof starts 
by expanding the microscopic evolution in a Dyson (or Duhamel) expansion. This expansion is standard and 
has been used many times in a related context. What might however seem odd at a first glance, is our complic- 
ated presentation of the perturbation series, involving tensor products of copies of the space £i = ,^(^i(J^)) 
(already introduced in Section [2. 1.3b . We use this formalism since we hope it makes the important estimates 
more natural and transparent. We first derive the expansion for Ga, G\, see the expressions l l3.17i[3A8l l. This 
part consists of purely algebraic manipulations (strictly speaking their validity is only established in Section 
I3.3.4l where we show that some series is absolutely summable). 

The analysis part comes in Section [Xil where we control ||G^||o- For reasons of readability, we first re- 
strict ourselves to 'bulk' polymers. Afterwards, the analysis is repeated with minor adjustments to include the 
boundary polymers; this is done in Section [3l6ll3.7l 

As mentioned above, we also need to establish the exponential ergodicity of the operator T, which is done 
in Section 13.51 The crux of this argument is to relate T to the Markovian approximation, which was already 
discussed in the introductory Section [L6l 

3.1 Dyson series 

Recall and rewrite the reduced dynamics Qt, introduced in Section [2. 1.21 
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By the Duhamel formula l|2.14|l in Lemma 123] we get 



men n^t.^ ^■ 



dh... dt2m Trp [Li,,(t2m) • ■ • ^1,,. (i2)ii,« (^1 ) (ps Pii)] (3.2) 



0<tl<...<t2m<t 



where, the RHS is xmderstood to be 1 for m = 0. Note that the RHS contains only terms with an even number 
of operators Li.k(s) because TrF(<I'Po) — whenever $ is a monomial of odd degree in creation/ annihilation 
operators. We use the Wick theorem to evaluate the expression in l|3.2l l. To this purpose, we use the formalism 
developed in Section l2.1.3[ as anticipated in Remark 12.21 For any ti,i — 1, . . . , 2m, we define a copy of the 
space £i (we do not aim to define "continuous" tensor products; all our formulas contain a finite number of t/s). 
In particular, we write 



5s...«)sii,K(i2)«)sii,«(ii)] (3.3) 
i™i^ti and it contracts the operators 



where, on the RHS, the operator T acts on operators in ■^{ti,....t2m} — '8 
such that those in .^f^ are on the right, then those on ^f^, etc., as in (|2.25t 

Let {u, v} be a pair of times with the convention that u < v. Then we define 

which can be written more explicitly: 

Ku^^ = -X^e^hiv - u) [n{D{v)) (g) C{D{u))] - X^e^hiu - v) \C{D{v)) ® TliDiu))] 
+ \^h{v - u) [C{D{v)) ® C{D{u))] + X^h{u - v) [n{D{v)) ® n{D{u))] 



(3.4) 



(3.5) 



where h{—s) = h{s), defined in Assumption lAl and D{s) ~ c^^^^Dc^'"^^. For later use, we note that ||i^u,ti||o < 
\'^C\h{v — u)\ for sufficiently small k (for example, C = 4\\D\\c'^). We view the operator K^^v as acting on the 
copies ® (that is, we should in fact write I ^u.v}[Ku,v] but, as for the operators G^, we prefer to drop the 
embedding operators /[•].) Let us now apply the Wick theorem and expand any contribution to the integral in 
l|3.2|l in contractions of the creation and annihilation operators. This yields 



rE[LL„(i2m))<8s 



E r 

7rePair(ti,...,t2™) 



(3.6) 



where the sum on the RHS runs over pairings tt, i.e. partitions of the times <i, . . . , t2m in m pairs {u, v} with 
the notational convention that u < v. This formula relies crucially on the fact that the interaction is linear in 
creation/ annihilation operators (nonetheless, it is easy to extend the proof so as to cover an additional small 
quadratic interaction). By plugging | |3.6|I into l|3.2|l , we obtain 



E 



0<tl<...<t2m<t 



dti . .. dt2m E ^ 

irePair(ti,...,t2T„) 



(3.7) 



where, for m = 0, the RHS is understood as 1. 



3.2 The evolution as an integral over time-pairs 

We will now rewrite Il3.7b in a convenient way. The integral over ordered <i, . . . , t2m, together with the sum over 
pairings, tt, on the set of times, is represented as an integral/ sum over unordered pairs {ui, Vi} with Ui,Vi G 
and i — 1, . . . ,m, such that 

Ui < Vi, Ui < . . . < Um (3.8) 

This is done as follows. For any pair {tr, tg} G tt with r < s, hence tr < tg, we let Ui — tr, Vi = tg where the 
index z = 1, . . . , m is chosen such that the Ui are ordered ui < U2 . • • < Um- We write w — {u,v} for a pair (with 
the convention that u < v) and w for a finite, possibly empty, collection of pairs. In the formulas below, we 
treat u, v as being implicitly defined by w. Given a Borel set J c R^, let E™ be the set of collections of m pairs 
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{ui, ?^i}i=i,...,m with Ui, Vi ^ J and satisfying the conditions I l3.8b . Let fimidw) be the Lebesgue measure on the 
corresponding subset - a simplex - of J^"' . We can then rewrite 



/ 



dti... dt2m ^ T 

0<ti<...<t2^<t 7rePair(ii,...,i2m) 



K 



{dw)r [(E) 



(3.9) 



where we abbreviated Ky^ := if„ „ for w = {u, v}. Let Sj denote the disjoint union 



m— J ' 



i.e. w e S,7 is given by w = (m, w™), m S N, € E™ with the convention Sj = {0} and we write |w| := m. 
Thus measurable functions F on Sj are collections {i^m}meN of measurable F,„ on E™. We let /i(dw) be the 
measure on Sj given by 

f fl{dw)Fiw) -^Y. I t^r,^{dw)Fra{w) (3.10) 



where we set /io((0, 0)) = 1- Note that the elements of Sj are naturally interpreted as sets, i.e. we write w E w 
to mean w e w™ for m = \w\ and w = w' U w" for the element in S'j with m = U with w € w whenever 
w £ w' or w £ w" . With these conventions the Dyson expansion (13.71 1 becomes 



(3.11) 



^[0,t] 



where, for w = the integrand is defined to be 1. 



3.3 Correlation functions and the Dyson series 
3.3.1 The contraction operator Ta 

To each macroscopic time t E In, we now associate a domain of microscopic times, 

Dom(T) = [\-^{t-1),\-^t] 
To a set A C /„ of macroscopic times, we then associate the domain 

Dom(A) = IJ Dom(T) 



(3.12) 



(3.13) 



TEA 



The contraction operator T[ ] defined in Section l3ll1 contracts operators so as to produce an operator in We 
now define a contraction operator Ta that produces operators in ^a- Let us first consider a finite family of 
operators Vt^ E where the indexed times ti satisfy ti < t^+i and ti E Dom(j4). Then we set 



Ta 





:= (E) It 


T 


i 


reA 





Vu 



j:tj GDom(r) 



(3.14) 



and we extend by linearity to the whole of lEi^ti, obtaining Ta ■ 'Ei^ti '-^ ^a- Iri words, Ta puts each operator 
into the right 'macroscopic' time-copy and contracts the operators within each macroscopic time-copy. 



3.3.2 The graph g{w) 

A set of pairs w E SDom(A) determines a graph Q {w) on In by the following prescription: the vertices t < t' are 
connected by an edge iff. there is at least one pair w = {u.v} inw such that 

u E Dom(T) and v E Dom(T') (3.15) 

We write Supp(5(w)) for the set of non-isolated vertices of G (w)/ i-e. the vertices that belong to at least one edge. 
If w; e SQofji(A) than Supp(t?(w;)) is a subset of A. In that case, we write Ga{ui) for the induced subgraph on A. 
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Let us denote the free S-evolutions 



=/,[e'(--i)^^], Ya= Yr. 

reA 



and = 7^[e-'^-^«], Ya = Yr 

TEA 



(3.16) 



We are now ready to state the connection between the Dyson expansion and the correlation functions Ga, G^: 
For Ac /„with|A| > 2, 



YaGaYa = 



n{dw) l[Supp(g(to))=A] Ta 



(3.17) 



^{dw) connected] ^4 



SDom(A) 

To check this, it is helpful to note first (with A = {ti, . . . ,Tm}) 



(3.18) 



YAE[Ur„ 



UrAYA 



^i{dw) Ta 



(3.19) 



which is a straightforward generalization of I I3.11II . Then, one observes that all pairs {u, v} with u,v £ Dom(T) 
make up It[T], and hence the operators Ga are integrals over all those w such that for each t E A, the graph 
Gini) has an edge {t, r'} containing t. This leads to (13.17b . Then, one verifies that G^ as given by Il3.18l l satisfies 
(|2.27l l. Since Il2.27b fixes G^ uniquely, (|3.18b is thereby proven. 



3.3.3 From pairings to C^: Pictorial representation 

We divide the time-axis into blocks (intervals of the form A~^[t — 1, r] with t e N)). We draw pairings as arcs 

connecting two times on the axis, either solid or dotted •, the reason to draw some of them differently 
is addressed below. The two times in a pairing, {u, v}, are denoted as Tegs'. This yields the upper picture in 
Figure 

Then, we single out blocks which do not contain a leg of a pair whose other leg lies in another block. Hence, 
such a block r can contain no pairing at all, or one pairing, or two pairings, etc... In the language of the graph 
G{w), this means that the time r has no edge to any other time; it is an isolated vertex. Such 'isolated blocks' 
are now, in the middle picture, drawn as -^^^ and we omit the pairs on them (imagining that they have been 

resummed and -^^^ stands for the sum). Blocks which are not isolated are drawn as • , but we still indicate 

the pairs on them. It is now clear that the pairings that were drawn with a dotted line in the above picture, are 
those inside an isolated vertex and only the other ones (drawn with a solid line) are reproduced in the middle 
picture. 

Finally, to arrive at the lower picture, we resum the pairings on the blocks • corresponding to connected 

components of the graph Q (w) . That, is whenever two blocks are connected by a pairing, we call them connected, 

and this induces a partition of the blocks • into connected components. Only this partition is indicated in 

the lower figure. Blocks belonging to the same (connected component) determine one polymer A. Alternatively, 
the operator is determined by summing all sets of pairs w that contribute to the connected component A. 



3.3.4 Bounds on correlation functions 

Up to now, we did nothing more on the Dyson expansions than straightforward algebraic manipulations. In 
what follows, we provide bounds. From the properties of the norm || ■ ||o discussed in l2.1.7l we get 

WC'aWo = \\YaG'XYa\U < I fi{dw) H \\K^\\^ (3.20) 

^Dom(A) 

To continue, we recall that Iji^tullo < \'^C\h(v — u) \ (Section l3.1| l and that ds{l + s)"\h{s)\ < oo by Assumption 
lAl We say that 'w spans A minimally' if the graph Ga{ui) is connected and no pair w can be omitted from w 
without losing this property. In particular, if w spans A minimally, than Qa (w) is a spanning tree. 
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Pairs determining for given A resummed 




r 



Figure 3: From the collection w of pairs {ui,Vi} to a collection A of polymers, A = {Ai, A2, A^, A4}. 



Lemma 3.1. Let \\h\\i := ds\h{s)\. For A c In, 



|G^||^<eCII''IUI^I /■ /i(du;)l[^ 



spans A minimally] 



n 



V — u) 



(3.21) 



Proof. We start with an appealing estimate that was the main motivation for encoding the pairings tt in the 
pair-sets w_. For integrable F, 

/i(du;) connected] I I < j Ai(dl«') spansAminimally] j /^(dw") l-F'lw' U w") | 

SDom(A) SDom(A) Soom(A) 

(3.22) 

To realize why this holds true, choose a spanning tree for the connected graph QAiw), pick a minimal subset 
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w' of the collection w such that GAiln') = 3' and use that 

Ai(d(w' U w")) = M(dw')M('^l«")- (3-23) 

We apply the inequality l|3.22|l to (|3.20l l with ll^wlU ^i^d we estimate the ^''-integral as fol- 

lows: 



V — u)\ 

.om(A) — 

m / Av,\^c\h(v,-u.-:)\ 



< 

^ (A^C|Dom(A)|||/.||i)" ^^A^c|DomU)|||/,.|h (3.24) 

~ ^-^ m! ~ 

meN 

and we conclude, since |Dom(A)| = A^^|A|. □ 
In what follows, we no longer trace explicitly the dependence on h and we simply estimate < C. 



3.4 Proof of bound (|2.51|) in Lemma 123 



In this section, we prove Il2.51b of Lemma 12.31 We start from Lemma 13.11 For each w_ that spans A minimally, 
Qa (w) is a spanning tree on A. Hence we can reorganize the bound (|3.21t by first integrating all w that determine 
the same spanning tree This amounts to integrate, for each edge of the tree, all pairs {u, v} that determine 
this edge. Hence we arrive at the bound 

||G^IU<e^l^l n ^(^'^') (3.25) 

trees spanning A {t,t'}gS{S) 

where is the set of edges of the tree ^ and, for r' > r, 

e(T,T') := j du J dv \^C\h{v - u)\, (3.26) 

A-2[r-l,r] A-2[r'-l,r'] 

and e(T', t) :— e(T, r'). Let At = t' — t and s = v — u, then 

(At+1)/A" 

e(T,T') < C y" ds\h{s)l At>1 (3.27) 

(Ar-1)/A2 
2/A^ 

e(T,T') < A^C / dss|;i(s)|, Ar = l (3.28) 







and, since ds{l + \s\)°'\h{s)\ < C , this implies the bound 

^ e„(T,T')<Ce (3.29) 

where e = |A|2niin("^i) (cfr. LemmaHJl and 

e„(r,r') := (1 + |r' - rj)" e(r, r'). (3.30) 

Starting from (13.25b , we bound 

A:A3to A:A3to trees .5" spanning A {r,r'}ef(5') 

< Y iCe)-^'^-^^^ n ^"(^'^') (3-31) 

trees 5^:roeV(5^) {r,T'}e£(3^) 
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where £{^),V{£^) are the edge, resp. vertex set of the tree 3^. To obtain the first inequahty, we used that 
rijr T'}e£{S)(^ + I''"' ^ ''"!)" d{A)°' for any spanning tree on A. The sum over trees in Il3.31b is estimated with 
the help of Lemma lA.ll by choosing 

1) the polymers S as unordered pairs {r, r'} and a{S) ~ 1. 

2) the adjacency relation S ^ S' iff. S" n 5' 7^ 0. 

3) the polymer weights ^({t, r'}) — {Ce)^^ea{T,T') 

The boimd Il3.29b plays the role of the Kotecky-Preiss criterion I IA7b and Lemma lA.ll yields 

sup {Cer\^'^^^\ n e„(r,r')<l (3.32) 

trees 5^:roeV(5^) {r.r'}e£(5^) 

For example, replace the restriction tq G by the weaker restriction (o'(5^) ~ {rg, ri} for some arbitrary ri 

and use ||A31 . Upon combining this with (|3.311 , the bound l|2.5Hl of Lemma [2. 3 1 is proven. 

3.5 Analysis of T: dissipativity and weak coupling limit 

Recall that operator T, acting on .^1 ( J^), that was introduced in Section l2.1.2l 

Tps = Trpb-'^^'^-Cps ® Pn)] (3.33) 
We exhibit the dissipative properties of T and we establish Statement 1) of Lemma [2.31 

3.5.1 Construction of Lindblad generator M 

To ease the presentation, we introduce some shorthand notation. Recall the one-dimensional spectral projections 
Pe corresponding to (see the discussion preceding Assumption |B). The set of differences of eigenvalues will 
be called £ = {e = e — e', e, e' € cr(iJs)}- Note that it is the set of eigenvalues of the Liouvillian Lg. By our 
non-degeneracy assumption, we associate a unique pair (e, e'), e, e' e (T(iJs) to any e e 5, e 7^ by e = e — e'. 
We then abbreviate = PeDPe' and we note that D* = since D = D*. We define the operator M, acting 
on^i(J^) 

Mp = -i[/fLamb,p]+ Y H-e)(e''D,pD;-^{D;D,,p}\ (3.34) 
with h{e) as in (11.9b , {A, B} = AS + BA, and the energy-shift ('Lamb-shift') operator 



ffLamb= ^ Im / dse'^^/i(s) i^^i?,, (3.35) 

which is well-defined by Assumption lA) . 

Because of our strong non-degeneracy assumptions, we can greatly simplify the form of M. Let us choose a 
basis for ^(J^?g) consisting of the one-dimensional projectors Pe and the operators |t/'e)('/'e' |/ with e 7^ e', where 
ipeecr{Hs) sre eigenfunctions of Hg. Then, inspection of ||3.34|| and the definition of the jump rates •) in lll.9b 
yields the following claims: 

Off-diagonal elements The operators \tpe){4'e' \ are eigenvectors of M with corresponding eigenvalue 

- i(^e, i?Lamb^e) " KA' , Hi^.mh^e') - (1/2) ^ {j [c' , c") + j{e, c")) (3.36) 

e" 

Assumption[B]iniplies that the second term is bounded away from and hence one has 

sup||e**^(|^e)(^e'|)ll <e-=* 

with c > . This corresponds physically to decoherence. 

Diagonal elements The space spanned by the projectors Pe is mapped into itself and we can identify the action of 



26 



M on this space with that of a ds X c^s-matrix (withds = dim J^) by setting A/(X]e Ai(e)^e) = J2eii-^^^)i^))Pe 
for /X e C^^. Then 

iMfi){e) = ^(e«j(e', e)M(e') - e'Ue)) (3.37) 

e' 

It is clear that, for k = 0, is the forward generator of a Markov process with state space o'(iJs)- We determine 
the spectrum of A^. Recall from the discussion following AssumptionlAlthat j{e', e) = whenever e > e', hence 

is a triangle matrix and its singular values are the diagonal elements. By inspection of | |3.37|I and Assumption 
|Bl it follows that all singular values lie in the region {z : Ke z < —Qm} with = ^^i^e-.ejtea J2e' ■? (^' save 
for a simple eigenvalue 0, corresponding to the eigenvector /i(e) = 6e,eo- By the spectral mapping theorem for 
semigroups with bounded generator, we then get that (T(e*^) C {0} U {z : \z\ < e^*^^ } 

Note that the right eigenvector does not depend on k, but the left eigenvector does. The Markov process at 
K = is called absorbing, with absorbing state eo. Physically, this property is due to the fact that the field F is in 
the vacuum state and cannot excite the atom S. We extract from this discussion what we need later on: 

Lemma 3.2. The ground-state projection Peo is an eigenvector ofe^^^ with simple eigenvalue 1. The rest of the spectrum 
gj^QtM i^^g 1^ jjjg region {z : \z\ < e~^^^'}for some gu > which does not depend on k. 

3.5.2 The Lindblad generator M and the microscopic model 

The relation between the the semigroup generated by the Lindblad generator AI and the microscopic model is a 
classical result in mathematical physics. It was discussed already in Section [L6l 

Proposition 3.3. For any i > 0, 

IIQ^ _ e-itis+A^tM|| < ^gCA^t|^|2a. (333) 

where a* — min(l, a) for a ^ 1, and |Ap"* = log |A||/or a= I. Using that Ls commutes with M, we get 

lime'^"''^sQ;,-2t = e'^^ (3.39) 

for any t > 0. 

For completeness, we review the simple proof of this convergence in Appendix [B] We now turn to the 
Proof of Statement 1) of Lemma \2.3\ Note 
i.) Lg commutes with and Af . 
a.) e''^^ is an isometry on (Mg)- 
Hi.) Ls{Pe) = for any e e (t{Hs), in particular Cq. 

Therefore, the spectral analysis of M (Lemma 13.211 implies that c^''^ '^Ls+m ^^^^ ^ simple eigenvalue 1 and all 
other spectrum lies inside a circle of radius e~^^'' . On the other hand, recall that T = Q\-2 and hence by 
Proposition l3.3[ T— e~'^ l^+m vanishes as A — > 0. We apply spectral perturbation theory of isolated eigenvalues 
to conclude that, for sufficiently small |A|, T has an isolated eigenvalue e®^ with corresponding spectral 
projector R such that statement 1) of Lemma [2.31 holds. Since T conserves the trace for k = 0, it follows that 
6't(A, k) indeed vanishes for k = 0; as A 0, it reduces to the eigenvalue of M, hence to 0, as well. Since 
T preserves positivity and the trace for k = 0, we can choose = 1 and 77 a density matrix, where ?/, fj were 
defined as the right, resp. left eigenvectors of T (see beginning of Section |Z3|| . For k e M, the operator T remains 
positivity-preserving and therefore 77, 77 can be chosen to be positive operators. □ 

For later use, we note that this perturbation argument also yields the bounds, for sufficiently small |A|, 

II77-III <0(|«;|), ||77-Peoll <qAp"*. (3.40) 



27 



3.6 Boundary polymers 
3.6.1 Dyson expansion 

We aim to write the analogue of the expansion (|3.111 for the expression of QnPs given in | |2.36|I . Recalling the 
definition | |2.10|I , we may write [/k in the same form as J2.14|l : 

/ dti...dt™ii,K(t™)...LL«(t2)ii,K(ii)ps®Po (3.41) 



'l<ti<...<tm<0 



with the definition 

= -£(1 ® $,(V'x , 0)) + 7^(l ® $_.(V'K , 0)), s e [-1, 0). 

In the same way we get 

U>,^J2 ("i)" / ■ • ■ Li,,{t„,) . . . ii,K (t2)iL«(ti) 

""^'^ n/A2<ti<...<t„<n/A2 + l 

where 

Li,Js) = -£(1 ® <i>(c'"l«l/^V>.,0)), s e (n/A2,n/A2 + 1)]. 
Thus we end up with a similar series as in (I3.2t and in (13.7b : 



(3.42) 



(3.43) 



(3.44) 



dii . . . dt2m Trp [^1,^(^2™) . ■ . Li^K{t2)Li,n{ti){ps ® ^h)] 

I<ti<...<t2™<ri/A2 + 1 

dti .. . At2m > ' T ® Ku,v 



-I<ti<...<t2™<n/A2 + 1 



7rePair(ti,...,t2m) 



{u,v}€:7r 



^(dw)T 



(3.45) 



-[-l.r./A2+i] 



where Ku,v is defined by the formula (|3.4| |, now with u, t; in [— 1, n/A^ + 1] instead of [0, n/\^]. Note that the 
choice of time-intervals [—1, 0] and [n/A^, n/A^ + 1] is somewhat arbitrary; we do this to display the boundary 
terms in a way that resembles the bulk terms as close as possible. 

3.6.2 Correlation functions 

First, we extend the definition of the graph Q (w) by allowing A c /„ and setting 

Dom(O) = [-1,0], Dom(n + 1) ^ [n/A2,n/A2 + 1] (3.46) 
We generalize Il3.16b to r = 0, n + 1 by putting 

yo = ?o = hm, Y^+i = J„+i[e'("/^')^^], ?„+i = /„+i[e-'("/^')^^] (3.47) 

and then defining Ya,Ya with A c /„ as in (|3.16b . With these definitions, the formulae (|3.17b , Il3.18b hold 
generally for A C In', 



YaGaYa = / A*(dw)l[SuppeA(H)=-4]'7A 



YaGaYa = / M(dl«) (to) connected] X4 



(3.48) 
(3.49) 
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3.7 Bounds on boundary polymers 

In this section, we prove the bound (|2.52b in Lemma [2.31 We first generalize Lemma [3. II to read 
Lemma 3.4. For A c /„ 



||G5 



Alio 



< Ce 



C\A\ 



f ^^{dw) I H \\K^\\A 1 



w spans A minimally] 



(3.50) 



Proof. To obtain bounds on the boundary pairings, we first note that (recall Dom(/„) = [0, n/\^]) 

' C\\h^(v)\ u e Dom(0),v e Dom(/„) 

C\Xhy, [n/X^ u e Dom(/„), w e Dom(n + 1) 

C|/iN(n/A2)| w e Dom(0),v e Doni(/„) 

CWij^f M,uGDom(0) 

C||i/'xi|P u, w e Dom(n + 1) 



We proceed as in the proof of Lemma [3. II and extract from each set of pairs w a minimally spanning subset 
w'. To bound the w" — w\ w' -integral in Il3.22b , we need an estimate on 



\Ku,v\\o < ' 



(3.51) 



Dom(A) UJGiO 



(3.52) 



Denote, for w = {u, v}, 



Xt<(w) 
Xixi(w) 



iMeDom(O) lu^Dom(n+l) 
iM^Dom(O) lueDom(n+l) 



IjieDom(O) li>eDom(n+l) 



(3.53) 

and Xt(liz) ■= n«,et« Xt(^) for l — k , x , N. For any w £ Doni(A), we can write w — w'^ U Wy^ U Wp^ such that 
w' G Dom(A n /„) and Xt(wJ = 1- Hence l|3.23b implies 

/ 



633 < 



\Sdo„ 



^l{dw) Y[ \\K, 



tullo I I Yl 

wew I \ te{x,xi,ixi} 



(3.54) 



The first integral is estimated in (13.241 1, and the three others in the same way using the bounds (I3.51I I. Thus, the 
i = K , X -factors are bounded by 



and the l =1X1 by exp(C|ft-ixi(n-/A^)|). Altogether we conclude that the product | |3.54|| , and hence l|3.52b , is bounded 

byCe^l-^l. □ 



To establish the bound I I2.52II , we closely follow the proof of Il2.51b given in Section 13341 First, we apply the 
definition of the edge factors e(T, r') — e{T' ,t) in II3.26II to the case where {t, t'} ^ /„; explicitly. 



e(T, t') 



C / ds|A/iK(s)| r = 0,T'e/„ 

Dom{r') 

C J ds\Xhy,{n/X'^ - s)\ T = ?i + l,T'e/„ 

Dom(r') 

C \h^{n/X^)\ T ^0,t' 



(3.55) 



29 



and then also ^^(t, t') := (1 + |r' — r|)"e(T, r'). These boundary edge factors satisfy 

e„(0,n+ 1) < (7, and ^ ea(T, r') < | AjC, forT = 0,?i + l (3.56) 

T'G/„ 

as follows immediately from the properties of/ix,/i,<i,/iix]in AssumptionlCl To get ||2.52| |, we first restrict ourselves 
to A with E A but rt + 1 ^ ^ (in particular, r = in ||2.52| |). Starting from Il3.50b and employing the definitions 
of Sair, t') above, we proceed as in (|3.31l l to derive 

1 ^ {Cl)-\^^^MAr\\G\\U < Yl (Ce)-l^(^)l n ^"(-'-') 

yle/„:A9O,n+10yl trees ,T:0eV(.T),n+l^V(5^) (T,T')e£-(,3^) 

where we used |An/„| > \S'{3^)\. For every tree ^ in this sum, we identify vertices ri < r2 < . . . < t,„ such that 
{0,Tj} e (S'(^) and rn > 1. Then the sum over 3^ is recast as a sum over m, choice of edges ti < t2 < . . . < t^, 
and trees growing out of them (including the possibility of no tree). This yields the bound 

m>l Ti<...<r,„ j = l y trees 5^:TjeV(5^)C/„ (r,r')e£(,S^) / 

^E E n2(C^£l^'e„(0,r,)< ^- 5^2(Ce)-ie„(0,r) <(7 (3.58) 

m>l Ti<...<T„, m>l ' \Te/„ / 

The first inequality follows by (|3.32b , since the trees are in the bulk, the last inequality uses H3.56b . The bound 
over A that do contain n + 1 but not and A that contain both n + 1 and 0, is analogous. Note however that for 
A that contain both n + 1 and we cannot extract an e factor for every edge since some trees contain the edge 
{0, n+1} (see in ll336l l). Therefore, the power of (e)"! in ll232l l is not \A\ - 1 but \A n /„|. 



4 Analysis of polymer model: Proof of main results 
4.1 Proof of Theorem ll.lt Approach to steady state 

To prove Theorem ll.il we first exhibit approach to a steady state for discrete times of the form t = nj}? and for 
a restricted class of initial states and observables (Section l4.1.1b . Then we eliminate these restrictions in Sections 
14.1.21 14.1.31 and 14.1.41 The reasoning in Section 14.1.11 is based heavily on cluster expansions and we advise the 
reader to read Appendix|A]before continuing. 

4.1.1 Approach to steady state for discrete times 

We start from the polymer representation (|2.71b with k = 0, hence = 0, 

E n«(^) (4.1) 

where Z„ = Zn{0, po, k = 0) depends on the initial state and observable via , fcx and the weights v{A) for 

An{o,n + i} 7^0. 

We write A ^ A' whenever dist(j4, A') < 1, and for a collection A, A ^ A' whenever there is at least one 
A E A such that A ^ A' (i.e., Suppyl ^ A'). These definitions coincide with those in Section lA.2l of Appendix lAl 
We separate each collection A into its boundary and bulk polymers by writing 



Zn^k^k^ E n"(^)pn,Supp.A (4-2) 

VA G .4 : A n {0, n + 1} ^ 
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where 

Zn,A'--= E n^(^) (4.3) 

A'*' A' 

Note that Zn,A' depends only on bulk polymer weights and that A consist maximally of two sets, one containing 
the element and one containing ri + 1. By the identity Il2.64b with k = 0, 

Z„j = Z„(l,?7(^Pn,K = 0) (4.4) 

By unitarity and the fact that Trs rj = 1, we have 

Z„,0 = Tr[e-'("/^')^(77 ® Po)e'("/^'^^] = Tr(7; ® Pn) = I (4.5) 

The quantity Zn^A' can be viewed as the partition function T„ of a polymer gas with polymer weights w{A) = 
v{A)l[AooA'] (i-e. it is of the form l|A5|l In Section rA.2b . For e small enough, the Kotecky-Preiss criterion l|A7|l is 
satisfied with S = Ce because of i2.75i with a = {a ^ will be used below), and hence Proposition lA.2l applies 
and yields 

log Zn^A'= E ^^W[A^A'] (4.6) 

where the truncated weights v"^ {■) , w'^ {■) are related to v{-),w{-), respectively, through the formula I IA61 in Sec- 
tion |A21 Comparing to the expansion of log Z„ j = 0, we get 

logZ„,^, = log^ = - E ^^i-^)MA^A'] (4.7) 

^"•0 AG'S,. 

By the bound | |A9|I in Proposition |A21 applied with Aq = A', we get immediately 

\ log ZnM <Ce\ A' \ (4.8) 

Next, we state 

Lemma 4.1. For \ X\ small enough and n — 0, the following limits exist 



lim V viA)Zn,A (4.9) 

1 — ^nn ' * 



n— >-oo 



- lim V viA)Zn,A (4.10) 



n— >oo 



Note that the number = (po) does not depend on O and z>c = z-^ (O) does not depend on po- 

Proof. To deal with the first limit, we recall the (partial) n- independence of polymer weights Il2.72|l , which 
enables, for n' > n, 

("') (") ("') (") y 

A A A:ma.^A>n A "-'^ 

in) 

where we abbreviated J2acI„-A30 n+ifA (i^ particular, the constraint e A is implicit in all terms) and 

A 

we indicated the n-dependence of polymer weights whenever its omission would cause an ambiguity. We use 
l|4.8|l and the bound 

(") 

J2 |f(A)|e""l'^l < m-"Ce, 

A:d(A)>m 
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which follows from G..76\ in Lemma IZ4l to conclude that the first term in (14. IH is bounded by Cm To bound 
the second term, we write, for max A <n, 

|logf^|< V lmaxsupp^>n|«^(^)| <Ce(l + n-max^)-" (4.12) 

The first inequality follows by inspecting l|4.7|l and the second follows by the bound 12.751 in Lemma [2. 4i which 
implies decay of cluster weights, as explained in Section lA.3l To conclude, we estimate the second term of (|4.1H l 
as 



(«) 



A ' A:maxA<ri/2 A:niaxA>n/2 



Using again Lemma IZH this is bounded by Cen^°'. 

The convergence in | |4.10|| is similar upon recalling the comparison l|2.73b between weights of polymers in /„ 
and /„' . □ 

In what follows, we now longer trace everywhere the factors e, e, thus making the bounds less sharp then 
possible. 

Lemma 4.2. For \ \\ small enough and n — 0, 

|Z„(0,po,'*-0)-fc^fc>,(l + z^)(l + z,,)| < Cn-" (4.14) 
and we identify fc^ (1 + ) = Tr p^. 

Proof. We start from ||4.2| | and we split the sum over A in five parts, of which the first corresponds to ^ = 0. To 
describe the other parts, let in general stand for subsets of /„ such that G A,n + 1 ^ A,\e\. Ay^ stand for 
subsets such that O^A, n + leA and, finally, write Ax , >] for subsets such that {0,n + 1} C A^x,. The splitting 
is 

{k^k^y^Zn = Z.nfi + ^v{Ay)Zn^A^+^v{Ay,)Zn^Ay,+ ^ w(Ak )w(Ax, )Z„^^^ UA>, + ^ w( A^ , xi ) , x, 

(4.15) 

We have already argued that the first term on the RHS is 1 and the second and third term converge to , Zx, 

respectively We split the fourth term on the RHS as Ea^.a>,:A^^a>, = Ea><.a>, ^ Ea><.a^:A^--a>, ■ Let us 
concentrate on the first term in this splitting, i.e. the unconstrained sum over , . Clearly, this expression 
should tend to Zx, . To see this we rewrite using (|4.7t 

v{A^)v{A^^)Zn^A^VjAy, - ^ v{AK)v{Ayi)Zn^A^Zn^Ay, = ^ u(At< )u(Axi )Z„_Ax -^n^A-, '7n(Ax , ) 

with 

g„(AK,A>,) =e~^-^e^""''(-^)^i^-*xi^i-^~-"»i - 1. (4.16) 
To bound g„ note first that (set |a|+ := max(a, 0)) 

dist(Ax,A>,) > \n-d{Ay)^d{A^)\+ (4.17) 

so that the bound IIA13b in Section lA31 together with Lemma lZ4l yields 

k^(-^)|l[>t~A,]l[A~A,] <C(l + |n-d(Ax)-d(A>,)| + )-"e|AJ. (4.18) 

Using |e'^^ - 1| < 7el"^l for < 7 < 1, a; G R we get from (|4l6t and (|418l 

\qn{Ay.A^)\ < (l + |n-d(Ax)-d(A>,)| + )-"e^'^l^xl. (4.19) 
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Furthermore, since for all , v4>c : 

(1 + |7i - d{A^) - d{Ay,)\ + )d{A^)d{Ay,) > cn, 

we conclude that 



(4.20) 



IRHSoflESH <C77-" \J2d{A^r\v{A^)\\Zn,Aje' 



C\iA^\ 



J2diA^rH^>^)\\Zn,. 



(4.21) 



The sums between brackets are estimated by Ce by first using l|4.81 and then Lemma [2. 41 
Hence, we have obtained the asymptotics of the first part of the fourth term in ||4.15|I : 



^ v{A^)v{Ay,)Zn,A^UA,, 



By similar, but rather simpler reasoning, we can now treat the other terms in | |4.15|| : 



E 



y^^v{A^)Zn^A^ 



v{AK)v{Ay,)Zn^A^UA, 



< Cn-°' 



< Cn 



'^v{A^^)Zn,Ay, Zy 
Ay, 



v[A^^y,)Zn.A^ 



(4.22) 



(4.23) 



(4.24) 



In fact, the upper line has already been obtained in the proof of Lemma [4. II This proves l|4.14b . 

Take now — 1, then Z„ ~ Trpo and, by the third remark following (|2.71b and using that 77 = 1 for 
K — Q, we have ky^ = \ and =0. Then, by the convergence established above, Z„ A:^ (1 + z^) and hence 
Trpo = fc><(l + ^k)- □ 

As promised, we comment on the case in which either or ky, vanishes, thus invalidating our expres- 
sion for the polymer weights given in (12.671 1. Let us assume for concreteness fcx = 0. Then all non-vanishing 
contributions to l l2.7Hl have S Supp^. All polymers weights v{A), A 3 are obtained by dividing a certain 
expression by fct< ■ Therefore, it is possible to redefine v{A) without the division and to omit k^ in the RHS of 
l|2.71b . It is straightforward to check that then the conclusions of the above lemmata still applies when fc^ (1 + Zk ) 
is replaced by z'^ where z'^ is defined by removing the division by fc^ in all its terms. The same remark applies 
to ky. 



4.1.2 Continuity of the polymer representation 

We envisage the situation that the model Hamiltonian and/ or initial state and observable depend on a parameter 
G 2? I, c C such that our assumptions are satisfied for any i' e 2?^. Note first that the numbers Znipo,0, k) 
are determined by a) the correlation functions /i, /ix , /ix , /im (we write h to denote any of the four correlation 
functions), b) the operators D, Hs,ps.o, Os in £M{Ji%) (we write Xs to denote any of these operators), and c) the 
parameter k. Indeed, the dependence on, for example, the Weyl operators, is hidden in the correlation functions. 
Below, we indicate the dependence on the parameter u by writing /j*^^-* and ^g"''- 

Lemma 4.3. Assume that 

1) The bounds of Assumptions^and^hold with the functions h{s) replaced by sup^g-p^ W'^H^)\- 

2) Assumption\B\is satisfied uniformly in v £ V^, more precisely, the gap gx and the constant Ct on the RHS of l|2.49b 
can be chosen uniformly in v. 

3) The operators X^'' are continuous in v and the functions hf"^ are pointwise continuous in v. 

Then, for sufficiently small |A|, the bounds on scalar polymer weights v{A) = v^^\A),A G /„ stated in Lemma \L^ 
hold uniformly in v and v'^^^ {A) is continuous in v,for any A e /„. Also the factors k^^ , fc^"-* are continuous in v. 
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Proof. By the dominated convergence theorem, one first estabhshes that the bounds on ||G5i||o hold uniformly 
in ly, and that the operators Tr , are continuous in ly (and the gap of T can be chosen imif orm in and then 
one passes from the operators Tr, G\ to the weights v^'^\ □ 

Note that we do not claim the continuity to be uniform in A as A — > 0. This can indeed not be deduced 
because of the rapidly oscillating factors e'^ in the definition of Ga, G\. 

4.1.3 From discrete to continuous time 

We extend the result of Lemma [4.2l to all times t, rather than times of the form n\^^ . Choose t > Q and imagine 
that we had started our analysis from Hamiltonian i^^H instead of H. This would modify the operators Xs and 
correlation functions h (notation as above) as 

D(^^=t-^D, H^^^^r^Hs, W\s)^h{r^s) (4.25) 

It is clear that all of our assumptions are satisfied for, say, I E [1,2]. Hence all conclusions of our analysis hold 
for this modified model as well. We focus on Z„ = Z„(0, po,K = 0) and let us indicate explicitly the dependence 
on £ by writing Zn \ By the result of Section 14.1.11 for sufficiently small |A|, the limit z!£ := lim„^oo Zn'' exists 
for any £ E [1,2]. Choose two values £1,^2 such that £1/(^2 G Q, then there is a subsequence of zi^^"* that coincides 
with a subsequence of zi^^-* and hence Z^''^ = Z^^K 

Since the correlation functions /i^^^ are Fourier transforms of L^-functions, they are pointwise continuous in 
£. Hence, Lemma [4.31 applies with v = £ and the polymer weights are continuous in £ e [1,2]. Since the bounds 
of Section l4.1.1l can be chosen imiform in £, we deduce that zlx^ is continuous in £, as well. Combining this with 
the arguments above, it follows that Z^ is in fact independent of £, and we can hence define 

(O)oo := (1/ Trpo) 1™ Tr[Oe~"^poe"^] (4.26) 

t—^OQ 

which enables us to proceed to the 

4.1.4 Proof of Theorem [Tl] 

The convergence l|4.26t is proven up to now for observables O and initial states po satisfying Assumption ICl To 
get point 1) of Theorem ll.il we show that the convergence holds for any initial density matrix po G ^i(^) 
and observable O e 2ITq. First, we note that vectors i/'x satisfying Assumption ICl are dense in L'^{M.'^). By the 
irreducibility of the Fock representation of the canonical commutation relations, this implies that vectors of the 
type ® ^(tAk)^ are dense in (see e.g. [25]). Hence, linear combinations of initial density matrices po for 
which (|4.261 holds, are dense in SSi (J^s)- On the other hand, by the unitarity of the dynamics, we have 

I Tr{Opt) - Tr(Op;)| < ||0||.g,(,^)||po - p^,||.^,(^) (4.27) 

and hence by a simple density argument the convergence (14.261 1 holds for any po G and O E Wa- 

Point 2) of Theorem ll.il Recall that here O is of the form (11.131 1 and ipy, satisfying 111. 15b . Without loss, we 
take IIOsll = 1- Then 

|(0)oo-Tr[(Peo'8f'o)0]| < \{0) 00 -k^\ + \k A -TT[{P,„®Pn)0]\ (4.28) 

< |fc>,z>,| + ||77-Peoll.^i(^s)l(f^,W(^>»)^^)l (4.29) 

< Ce\k-^\+C\X\^'^' (4.30) 

To obtain the second inequality, we used that (0)oo = fcxi (1 + Zy, ), as obtained in Section [4.1.1l and the definition 
(12.701 1. The third inequality was established in the proof of Lemma [4.11 (bound on Zy,) and in ( 13.4011 (bound on 

Point 3) of the Theorem was established in Section [4.1. II for discrete times and is immediately extended to 
continuous times by the reasoning in Section [4.1.3l Hence the full Theorem 1 1.1 1 is now proven. 
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4.2 Photon number bound 



To derive a photon number bound, we are interested in Tr[/3(e''^], that is, Z„ with = 1 and general po- 
However, it is more convenient to investigate first a sHghtly different quantity; namely 

p{k) := lim - logTr[(ry (g) e^^)p„/x2], Rt) = V ® Pn (4.31) 

with ?7, ?7 as in I l2.50b . We call this quantity a 'pressure', since it is the logarithm of Z„, which we view as a 
partition function. The advantage of the definition (14.311 1 is that we can establish right away 

Lemma 4.4. For \ k\ small enough, the pressure p{k) exists and is analytic in A 

Proof. Let Pn/x^ be as above, then by Il2.63b and Il2.64b , 

Tr[(77(ge'=^)p„/;,2] = Z„(,7®l,r7 0Po,«) = Trsi?Q„/A2i? = e"«- ^ H "(^) (^.32) 

Starting from Il4.32b and using the cluster expansion results in the Appendix section lA.41 we get thatp(K;) exists 
for e (hence A) sufficiently small. To verify that Section lA.4l is applicable, we need Lemma |Z4| and the translation 
invariance of polymer weights, see below ll2.71b . 

We argue that p{k) is analytic in n by proceeding as in the proof of Lemma |431 Note first that the parameter 
K enters only (at least for bulk polymers, which are the only ones concerning us here) by multiplying two of the 
four terms in Ku.v by c", see I l3.5b . Hence, by the Vitali convergence theorem, the operators GJ^, T, and also the 
scalar polymer weights v{A) are analytic in n. Since analytic polymer weights imply analyticity of the pressure 
by the Vitali convergence theorem (see Section lA.4.Hl the lemma follows. □ 

Of course, given sufficient infrared regularity, we expect p{n) = 0. This will be established in Section 14.2.21 
by introducing and removing an infrared cutoff. However, we would like to draw attention to the fact that there 
is in principle a more straightforward way. One could establish that lim„^oo ^ logTr[e'*^p„/;^2] is independent 
of Pa, and then use the existence of a ground state ^^gs with the property that (^gs, e^^^gs) < oo to argue that 
lim„^oc log Ti-[c'^^ pn/x2 ] must vanish in general, since it vanishes for po = 1 5'gs) (5'gs | . We choose not to exploit 
this approach so as to keep our presentation as self-contained as possible. 

Once we have p{k) = 0, we still have to exclude that Tr[e'^^ pt] grows subexponentially in t, and this is done 
in Section l4T3l 



4.2.1 Infrared cutoff and spectral translations 

In this section, we introduce an irtfrared cutoff 7 and we argue that the pressure p{k, 7) is continuous in 7. 

We implement the cutoff by translating the form factor in the spectral parameter of the one-photon Hamilto- 
nian, i.e. the multiplication with \q\ on L^(R, dq). We write q = ujq with lo = q and q G S"*"^, the sphere equipped 
with the surface measure dq. Consider the isomorphism of Hilbert spaces 



J : L^(M^dg) ^ L\M.+ ,L\E>''~\dq),duj), 



(4.33) 



such that {Jij:){u]) is the function g ^ cj 2 ip{ujq) in L^(S ^, dq). We define the '7-translated' form factor 
by specifying Ji/)''''^: 



(J0(^))(c.) := 



J(j)(u) — 7) w > 7 



w < 7 

In general, the correlation function h{-) can be written as 



for7 > 



h{t) = / dqe-'l*l*|<^(g)p = / dc.e--*|| J0(c.)L2(s.-i) 



(4.34) 



(4.35) 



It is immediate that hf^^ {t) = e '''*/i(t) and hence in particular 

/i(7i)(i) _ /j(72)(^) < i-e't-^i--^^)* \h{t)\ 



(4.36) 
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We want to apply Lemma 14.31 to conclude that the polymer weig hts are continuous in 7, and con- 

sequently (see Section lA.4.1l l also 7 i-> p{k, 7) is continuous. In fact, the conditions of Lemma [4.31 demand that 
we also check continuity of the functions h^,hy,, /im, but this is not necessary as long as we only use bulk poly- 
mers (since those do not depend on , hy^ , h^). Indeed, the pressure p{k,'j) above depends only on the bulk 
polymers and hence we conclude that p{k, 7) is continuous in 7 e [0,1]. 



4.2.2 Pressure vanishes 

The advantage of having a sharp infrared cutoff is that we can bound the expectation of the number operator in 
the state in terms of the expectation in the state ^'q. Indeed, let us decompose the number operator N as 

N^N^'' + N<'', N^^ -.^ [ dqa^aq, N<^ := [ dqa^aq (4.37) 

J|g|>7 "'|9|<7 

Thenbecause of the cutoff, [iV<''^,i?] = 0,andwecanboimdiV^'^ < 7" liJp- Choosing I* such that (^'o,iV5'o) < 
C and {^q,H'9q) < C and H^oH = 1, we estimate 

7(*t,^*t) < 7{-9o,N<'<^>o) + {^t,HF^t) (4.38) 

To get an estimate on {^'t^Hp'^t), we write Hp = H - Hs - XHi, we use that (^'t,iJ*() = (*o,-H^*o>, the 
boundedness of Hg and the estimate l|1.7|l to obtain 

(*t, ilF*t) <C+ |A|C"(*t, Hp^t) (*t, i^F*t) < C" (4.39) 

for A sufficiently small. By using the spectral theorem and Jensen's inequality, we can now bound (*t,e'=^^'t) 
with K < from below and obtain 

lim -log(5' - ,6'^^^' " ) = 0, K<0 (4.40) 

n— foo 71 a2 a2 

We want to deduce that also p as defined in (I4.31t vanishes. This indeed follows immediately. Firstly, since the 
operator 77 in the definition of p{k, 7) satisfies — =C'(|k|) and for real k, it is a positive-definite operator, 
we can find constants c, C such that cl < fj < CI for sufficiently small, real k. Secondly, the density matrix 
Pa = 77 (g) is a finite sum of terms j^'o) (^'ol with satisfying all conditions necessary for the estimate l|4.38|l . 
Hence we have shown 

p(k,7) = 0, k<0,7>0. (4.41) 

However, by analyticity in k (Lemma l4.4b . this holds for sufficiently small |k| and by continuity of the pressure 
in 7 (established in Section l4.2.1b , it holds for 7 0, as well. 



4.2.3 Proof of Theorem 11.31 Finite-size corrections 

The fact that p ~ tells us merely that the number of emitted photons grows slower than linearly in time, a 
result that has been established with other techniques as well. We aim to prove that it does not grow at all. Let 
us denote the quantitiy in I I4.32II by ^„_0, in analogy to the definition in Section |4. 1 . 1 1 (the only difference is that 
now K ^ 0). We invoke the discussion in Section lA.4i in particular the formula l|A17b which, applied to the case 
at hand, reads 

logZ„,0 = - Yl min(d(^)-l,7i)t;^(^) (4.42) 

^eSoo :min Suppy^l=l 

where we also used that p — 0. We bound J4.42|l as 

|logZ„,0| < C'7ii-™("'i) J2 rf(yt)"|f^(y^)| < Cni-™("'i) (4.43) 

yle S oc :min Supp.A= 1 

by Proposition lA.2l Next, we show that this bound remains valid for arbitrary initial state satisfying Assumption 
ICl We proceed as in Section l4.1.1l obtaining 

^" fcKfc. E fn«(^)|%^ (4.44) 



A€<Bl \AeA 
VAe A: An{0,n + 1} ^<ll 
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In contrast to the situation in Section l4.1.1[ Zn 7^ 1 since 7^ 0. However, the second equality in <|4.7|) still applies 

z ~ ^ 
and, in analogy to l|4.8|l we get | log "^^"^^"^ | < eC|Suppyi| and then also 

1^1 <C, (4.45) 

by combining (|4.44b and Il2.76b and making e sufficiently small. 
Hence we have in general 

l^n|-|#^ll^„.0l<C'e^""™'°"' (4.46) 

Finally, we mimic the reasoning in Section l4.1.3l to extend the bound (14.461 1 to continuous time. This finishes 
the proof of Theorem [O] 



A Cluster Expansions 

We present some standard results on cluster expansions. We opted not to present a very general exposition, 
but instead a framework that is maximally adapted to our problem. A slight exception is Section [A. II where 
we state a combinatorial lemma that is used in several places in our paper. In that section, the 'polymers' can 
be anything, but in the remaining sections of this appendix, we assume them to be subsets A of N, just as in 
our construction, and although we call the polymer weights w{A), the reader can safely confuse them with our 
concrete polymer weights v{A) defined in Section l23l In the sections other than Section lA.il we also freely use 
the notation introduced in Section [2.1. 51 This appendix follows very closely the presentation in l30l , which uses 
the convergence criterion by I22J . 



A.l A combinatorial lemma 

Let .5^ be a countable set. We call its elements S £ ,y polymers. A function w : y ^ C will be called a polymer 
weight. Furthermore, let ^ be an adjacency relation on y i.e. ^ is a symmetric and irreflexive relation. Given a 
finite collection of polymers S C y we define 

^(^5)^= E n Ms.^s,] (Ai) 

where the sum runs over ©^(5), the set of connected graphs with vertex set S, and fS'(^) is the edge set of the 
graph Finally, for a collection S and a polymer 5*0, we write S ^ Sq whenever there is a S* G 5 such that 
S ^ Sq. Then, we have the following combinatorial result 

Lemma A.l. Assume that there is a positive function S a{S) such that the so-called "Kotecky-Preiss criterion" 

^ \w{S)\e<'^ <a{S') (A2) 

S:S~S' 

holds. Then 

J2 HS)l[\w{S)\<a{So) (A3) 

and 

J2kiSU{So})ll\wiS)\<e^^'''^ (A4) 
5 ses 

where the term corresponding to S = % is understood to equal 1. 

An easy way to prove (|A4b is via induction in the number of elements of the collections S. The claim IIASI l 
then follows by choosing Si £ S such that Sq ~ 5*1 and applying (|A4b with Si in the role of 5*0. 
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A.2 Logarithm of the partition function 

We make the setup introduced in Section lA.ll more concrete 

1) The polymers are now subsets of /„ = {1, . . . , n}. They are denoted hy A, A' , . . .. 

2) The adjacency relation is defined tohe A ^ A' 4^ dist(A, A') < 1 with dist(A, A') = inf reA.r'eA' k - t'|. 

3) Collections of polymers are denoted by A. We use the sets of collections *B„,*B^ introduced in Section 

ma 

We start from a polymer-gas representation of some partition function T„: 

AG'Sl^AeA Ae<S„AeA {A,,Aj}cA 

where the last product runs over pairs in the collection A, and the summand is defined to be 1 if ^ = 0. 
We define the truncated weights oi A<E 25„ as follows 



T 

W 



{A) := E (-1)'^^^^' n n -(^) (A6) 

'Se&^(A) {A,,Aj}e<f{'^) A,eA 



where the sum runs over <3''{A), the set of connected graphs with vertex set A, and is the edge set of the 
graph We call A a cluster whenever the graph on A with edge set {{Ai, Aj}, Ai ^ Aj} is connected. Hence, 
if A is not a cluster, then w'^{A) = 0. 

Next, we state the basic result of cluster expansions (cfr. (eq. 4) in II30I ). 

Proposition A.2. Assume there is S > such that the 'Kotecky-Preiss criterion' 

e^\^\w{A) < S\A'\ (A7) 



AcI„:Ar^A' 



holds. Then T„ ^ and 
and, for any A^ c /„ 



Y \w^{A)\<6\Ao\ (A9) 

A'^'Sn-Ar^Ao 



This proposition is a direct consequence of Lemma lA.ll It is clear how it applies to the present paper. The 
polymer weights w{A) are v{A) as constructed in Section |Z31 and the Kotecky-Preiss criterion is (12.75b where 5 
should be chosen not smaller than Ce and not larger than c (C, c as in Il2.75l l). 

A.3 Decay of cluster weights 

We assume now that the Kotecky-Preiss criterion is satisfied in a stronger sense, namely: 

Y e'\''\w^{A)<5\A'l (AlO) 

where 



and we show that this yields some decay in the cluster weights. Proposition |A2] gives 

E ha (-4)1= E Yld{Arw'^{A)<5\A,\ (All) 

Ae'&^-.A^Ao Ae'&^-.A^Aa AeA 

If is a cluster then X^ag^a ^(^) — ^^(-4)' ^^'l if additionally all polymers A with nonzero weight have d{A) > 1 
(this is our case), then 

W d{A)" > d{A)" (A12) 
AeA 
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Hence, we use <Allb to derive 

UA)>rn\w^{A)\ < CS\Ao\{l + m)-'' (A13) 

A.4 Pressure 

In the previous sections, n was a fixed parameter. Now, we consider different values of n. If A c /„ C /„' 
(n' > n), we assume that w{A) does not depend on whether we view A as polymer in /„ or /„/ . Moreover, we 
assume that the weights are translation-invariant in the sense that 

w{A) = w{A + a) (A14) 

where A + a := {t,t — a £ A}. An immediate consequence is that the pressure 

p :— lim n~^logT„ (A15) 

exists and is given by 

p= ^ w^iA) (A16) 

.Ae SB oo , min Supp>t= 1 

where 03 oo = U„*B„ is the set of finite collections of finite subsets of N. Indeed, the absolute summability of the 
RHS of IIA16b follows from Proposition [AjJl whereas (we set |a|+ := max(a, 0)) 

p-n^MogT„ = ^ 

.Ae SB oo :min Supp>t= 1 

which vanishes as n — > oo be the dominated convergence theorem, since the factor between brackets converges 
to for any value of d{A). 



n-\n+l- d{A)\_ 



w^{A) 



(A17) 



A.4.1 Continuity and analyticity of the pressure 

Often, one encounters the situation where the weights w{A) depends on a parameter ly G Vi, c C, w{A) — 
w'^"^ {A) such that the Kotecky-Preiss criterion ||A7| | is satisfied even when w{A) is replaced by sup^gp jw^'^' (^)l- 
Then 

1) If the weights w^^^ {A) are continuous in e Vy, for any A C N, then the pressure is continuous in e Vi,. 

2) If the weights w^^^ {A) are analytic in i/ e Vjj, for any A c N, then the pressure is analytic in e V^. 
Starting from the absolute summability of the RHS of (|A161 , the first claim follows from the dominated conver- 
gence theorem, and the second from the Vitali convergence theorem. 



B Van Hove limit 

In Section l33l we introduced the Lindblad generator M and we stated and exploited the claim (Proposition l3.3t 
that 

WQt - e-"^^+^'"^|| < c^^'*|A|2"* (Bl) 

Below we derive this claim. The derivation is independent of the rest of the paper, except for the representation 
of the effective dynamics Qt- 

Qt^ f M(dw)r| <E) kA (B2) 

and the bound ll-ftTn.tjIU < A^|/tft) — u]\. Moreover, we freely use AssumptionlAl i.e. dt(l + < C. We 

set K = in this Appendix. Including a nonzero value for k does not affect the reasoning. 
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B.l Alternative construction of M 

We define the operator e ^ (J^s)) using the notation of Section|3l 

A2f,_, e-'''^^r[i^„,.]e'"^^ (B3) 

For the sake of clarity, we write out Fg explicitly (recall that k = here) 

Fs = -h{s)'R{D)c-'''^^ C{D) - h{-s)£{D)c-''^<'TZ{D) 

+ h{s)£{D)c-''^''£{D) + h{-s)'R{D)c-'"'^^'R{D) (B4) 

Introduce spectral projections Pees of the Liouvillian Ls- The range of Pe=o has dimension dim ,3% and the Pg^o 
are one-dimensional. The connection between the operators Fg and the generator M given in Section 1331 is 

/•OO 

M = ^ / dse'^^P^FsP^ (B5) 



This can be checked by direct computation. The four terms in llB4t are matched with the expression (|3.341 
as follows: the first two terms give the first term imder the integral (with k). The two other terms yield the 
commutator with Hiamb and the second term under the integral. Note also that the part of M corresponding 
to the e = term on the RHS of dBSb is isomorphic to the operator M discussed in Section |33]lJ and the parts 
corresponding to e 7^ were discussed imder the heading 'off -diagonal elements' in that same section. 

B.2 Emergence from the Dyson expansion 

We start from the Dyson series (|B2) and we split the set of collections of pairs ^[o,t] iuto leading Sj^ and 
non-leading EjjJ subsets, as follows 

^o,t] ■= e S[o,t], Vi = 1, . . . , - 1 : < u^+i}, SjjJ := T,[o,t] \ ^o,t] (B6) 
(note that e S{g ^j). Define 

Qt,i ■= J fJ.{dw)T {(^wew^w} 

00 „ 
m=0 

Its Laplace transform is 



0<...<Mi<Di<«i + l<...<t 



dt e-''Qt,i =(z + iLs - A^F,) ' (B7) 



+ 



where Fz = /^^ dt e *^Ft. By inverting the Laplace transform we get 
Lemma B.l. 

where a» = min(l, a) for a ^ 1, and |Ap"' — A^| log |A||/or a = 1. 
Proof. We start by writing 

Qt,i - e-'*^^+^'"^ - ^ dz{Qzj - g,,o)e*^ (B9) 
where C is a contour in the complex plane, of the form C = iR + AX^ with A large enough and 

Qz,o = {z + iLs-XHl)-\ 
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Indeed, since the spectrum 8 of is is real and ||-Fz|| is uniformly bounded for Rex > 0, the function Q^^i is 
analytic in the region Re z > if A is large enough and so (IB9b holds. 
Next, simple algebra yields 

Qz,y - Qz.o = X'^QzAF, - M)Q,,o (BIO) 

and from dBSII 

M = Y,PeF-uPe- (Bll) 

We decompose = Z^ees /c where {z e C, |z + ie| = niine'gg |z + ie'|}. For z e we write 

- M = {f^,, - Af) + (f, - F_i,) K + W,. (B12) 

The lemma follows if we show 

A2 V / dz||g,,i(K + W^s)Q..o|| < CA^"-. (B13) 

We need the bounds, for z G C, 

llQ.jll < C(dist(z,-if))-i (B14) 
IIQ.j-Psll < Clz + iej-i (B15) 
for i ~ 1, 0. The bound for j = 1 follow by straightforward estimates on the Neumann series 

— 1 n 

(z + {L^-\^F,^ = (z + iLs)"'^ (A'-F.(z + iLs)"')" (B16) 

n=0 

using in particular \\{z + iLs)^^|| < C(dist(z, — i£))~^ < C|A|^^. The j = case follows similarly. Using jBlH 
we have 

K = Pe'F-iePe" 
{e',e"me,e) 

so combining with (|B14b and IIB151 

||0..iKQ.,o|| <C V — — , forzeCe. 

The K contribution to l|B13t is then bounded by 

CA^ V / dz , , . ,1 , . „ <CA^|log|A||. (B17) 

As for the we have 



Z + 1£ Z + 1£' 



\\Fz~F-^e\\< Ath{t)\e-~^^ -e'^^\. (B18) 
Jo 

Using dt + t)" < C, we obtain for Re z > 

\\Fz - F-»e|| < Cmin{|z + ier'"("'i), 1}. (B19) 
Hence the contribution to IIB13b is bounded by 



(B20) 



Y f dz ^— min(l, |z + i£|™"("'i)) < CjAP"- . 

(lBT3l follows now from l|BT7l l and l lB20ll . □ 
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It remains to estimate the contribution of the non-leading pairs: 
Lemma B.2. For any i > 0; 



sup / dii;||r{®^eH^^} II <e*<='|Ap'™"(i.«) (B21) 



Proof. We start with a bound analogous to I I3.20II : 

j dw\\T{<)^n,en,K^]\\< j Aw\{\^C\h{v~u)\ (B22) 



By definition, every w e SJjJ has to contain at least two pairs (u, w), (m', w') such that u < u' < v. Choose the 
first two such pairs and integrate over the coordinates of all other pairs, proceeding as in l|3.24t . This yields 

J dwll CX''\h{v -u)\< Ce^'^*ll''lliA4g(t) (B23) 



with 



q{t)^ J dudvdu'dv' \ h{v - u)h{v' - u')\. (B24) 

0<u<u' <v<t,u' <v' <t 

q{t) is estimated by first performing the integral over v', which gives a factor \\h\\i, and then the one over u', 
giving |v — u|: 

q{t) < \\h\\i J dudu I - It) 111) - m| 

0<Ji<u<t 

< ||/i||i^'dz.(^^max Jl + |.|)i-")^°°d,s(l + |.^ 

< ||;i||iC<2-2min(l,a) (g25) 

Taking < t < |A|^^t, this yields the lemma. □ 
The boimd jBH (and hence Proposition |3]3]l follows by combining Lemmata lB.2l and lB.ll 

C Infrared Regularity 

We clarify how our Assumption |Al formulated in the time-domain, relates to infrared regularity of the form 
factor ip. First, we state the auxiliary 

Lemma C.l. Let f : (0, oo) ~^Rbe a measurable function of compact support and such that for some < 7 < 1, 

|9'7MI <Cou;'^"i-", n = 0,l (CI) 

Then, for any /? < 7, 



/(w)e'*"dw 







Proof. Define, for < /3 < 1, 



1/1^:=/ dc. sup e-^|/(c. + e)-/(c.)| <oo (C2) 

Jo 0<e<l 
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Then, we claim that 



/(a;)e*'"da; 



<j^^—^{\m\-^ + C^-'\tr) (C3) 



Lideed, writing (say t > 0, the case t < follows by replacing / by /) 







the claim follows. To verify that |/|/3 < oo for any /3 < j, we estimate 

e-^\f{oj + e)-f{u)\<{ 

\dfiu;')\ 

and then use l ICll to bound this by Cw^^+"''~^, which is integrable for any /3 < 7. □ 
Using the above Lemma, we can state a sufficient condition on the form factor (/> for Assumption lAl to hold. 

Lemma C.2. Fix 7 > (not necessarily 7 < 1) and let [7J he the highest integer not larger than 7. Let fq{u}) := 
a;'^~^|(/)(a;g)p, with q e S''^^ and assume that 



|a:/^MI<c, Vc.e(o,2]: la^v^MI < cc.^-", (C4) 

/or any n = 0, 1, . . . , [7J + 2 and uniformly in q. Then, Assumption\A\holds with a < 7. 
Proof. Let 6* : M+ ^ [0, 1] be a C°° function such that 

^h4; -"tl (C5) 
10 w > 2 

and write / = /^^^ + /^^^ with /^^^ = 6*/ and /^^^ = (1 — 0)/ where we omit the parameter q (all bounds are 
uniform in q). Then J dcj/^^)(w)e^'*" decays as t^^fi+^ by the first bound in ||C4|| and partial integration. To 
estimate / dw/'^^^(w)e^'*", we use the second bound in iC4i and partial integration [7J + 1 times to extract a 
factor Then we use Lemma fC.ll but with 7 — [7J < 1 playing the role of 7 in Lemma [Cli to extract the 

additional decay t^^ with /3 < 7 — [7J . □ 

Likewise, one can derive from Lemma [C . 1 1 sufficient conditions on 0, i/j^ , to satisfy AssumptionlCl 
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